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FOREWORD 


This  is  the  final  report  on  Contract  AF29(601)-2596,  ARF 
Project  K199.  "The  Investigation  of  Silo  and  Tunnel  Linings",  which 
commenced  in  February  I960.  The  work  was  administered  under  the 
Air  Force  Special  Weapons  Center,  with  Mr.  C.  Wiehle  as  project 
monitor 
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and  T.  A.  Zaker.  Particular  acknowledgement  is  made  to  M.  A.  Salmon 
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important  contributions  for  the  remainder  of  the  contract 
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ABSTRACT 

Some  of  the  aspects  of  the  problem  of  determining  the  response  of 
underground  cylindrical  tunnel  and  silo  linings  to  blast-induced  loading  are 
presented  in  this  report.  The  studies  include  the  effects  on  silo  linings  of 
differential  horizontal  soil  displacements,  vertical  skin-friction  forces 
generated  by  vertical  compression  of  the  surrounding  soil,  and  axial  response 
of  a  silo  with  integral  cover  to  the  blast-induced  cover  loading.  (U) 

A  theoretical  method  for  estimating  the  forces  acting  on  a  cylindrical 
structure  subjected  to  air  blast-induced  ground  shock  is  presented.  The 
radial  forces  on  the  surface  of  the  cylinder  consist  of  the  undisturbed  free- 
field  pressure  component  plus  pressures  proportional  to  the  radial  velocity 
and  displacement  of  the  structure  relative  to  the  free-field  particles.  The 
effect  of  virtual  mass,  foundation  modulus,  and  damping  on  the  response 
was  studied  by  obtaining  solutions  appropriate  for  tunnel  and  silo  linings. 
Static  tests  on  three  metal,  cylindrical  tunnel  models  buried  in  dense  Ottawa 
sand  are  presented  and  compared  with  the  theoretical  method.  (U) 
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INVESTIGATION  OF  SILO  AND  TUNNEL  LININGS 


I.  INTRODUCTION 

This  is  the  final  report  on  Air  Force  Contract  AF  29(  601)  -2596, 

"The  Investigation  of  Silo  and  Tunnel  Linings",  initiated  on  February  26,  i960 
and  extended  on  January  5,  1961.  This  report  covers  the  work  performed 
by  the  Armour  Research  Foundation  during  the  course  of  the  contract. 

To  prepare  adequate  economical  designs  for  silo  and  tunnel  linings 
to  resist  blast  induced  loads,  an  understanding  of  the  behavior  of  these 
structures  in  their  soil  environment  is  necessary.  The  objective  of  this 
program  was  to  develop  methods  for  the  analysis  of  the  response  of  silo  and 
tunnel  linings  to  blast  loads.  It  is  anticipated  that  the  results  of  this  study 
can  be  incorporated  in  a  rational  method  of  design  for  these  structures. 

The  air-blast  pressure  traveling  over  the  ground  surface  produces 
stresses  and  displacements  in  the  soil  in  both  vertical  and  horizontal  direc¬ 
tions.  Since  the  resistance  to  deformation  of  an  underground  structure 
differs  from  that  of  the  soil  it  replaces,  a  knowledge  of  the  state  of  stress  in 
the  free-field  is  insufficient  to  permit  direct  determination  of  the  loading 
on  the  structure.  That  is,  some  consideration  must  be  given  to  the  major 
problem  of  estimating  the  effects  on  loading  of  the  soil  -  structure  interaction. 

Consider  the  effects  of  air  blast  induced  ground  shock  on  a  cylinder 
of  soil  which  a  silo  or  tunnel  lining  replaces.  The  lateral  surface  of  the 
cylinder  will  be  subjected  to  certain  stresses  which  will  produce  deformation. 
These  stresses  and  deformations  are  known  to  the  extent  that  free-field 
effects  can  be  predicted.  In  general,  the  deformations  will  be  such  that 
originally  circular  cross  sections  of  the  cylinder  will  become  elliptical  in 
shape,  the  axis  of  the  cylinder  will  become  curved,  and  there  will  be  a 
shortening  of  the  cylinder  in  the  axial  direction.  The  lining  which  replaces  the 
soil  will  tend  to  deform  in  the  same  way.  The  elliptical  deformation  of 
originally  circular  cross  sections  will  produce  circumferential  bending  stresses 
while  the  uniform  component  of  the  radial  pressure  will  cause  compressive 
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circumferential  stress.  The  curvature  of  the  axis  of  the  lining  will  produce 
longitudinal  bending  stresses  in  the  lining  acting  as  a  beam  and,  since  the 
lining  will,  in  general,  be  more  rigid  in  axial  compression  than  the  column 
of  soil  it  replaces,  longitudinal  friction  forces  between  the  soil  and  the 
lining  will  produce  direct  axial  stress. 

To  simplify  the  analysis,  each  of  these  effects  is  considered 
separately,  as  follows. 

1.  An  analytical  method  for  estimating  the  interaction  effects  on 
the  response  of  a  circular  cross  section  of  such  cylindrical 
structures  as  silo  and  tunnel  linings  is  presented  in  Section  II.  In 
this  analysis,  an  attempt  is  made  to  account  for  the  soil- structure 
interaction  by  introducing  load  components  proportional  to  the 
displacement  and  velocity  of  corresponding  points  in  the  undisturbed 
free-field.  Inertial  effects  arising  from  the  disturbance  of  the 
free-field,  caused  by  the  presence  of  the  structure,  are  accounted 
for  by  the  addition  of  a  virtual  mass  to  the  actual  mass  of  the 
structure.  The  constants  of  proportionality  relating  relative 
displacement  and  velocity  to  load  are  estimated  on  the  basis  of  con¬ 
sideration  of  certain  elastic  media  properties.  The  effect  on  the 
dynamic  response  to  variations  in  the  constants  are  then  studied  to 
determine  the  sensitivity  of  the  results  to  the  parameters.  The 
analysis  should  be  of  value  in  estimating  the  relative  importance 
of  dynamic  effects  on  the  response  of  silo  and  tunnel  structures,  and 
in  establishing  the  range  of  applicability  of  simple  design  criteria 
based  on  static  loadings.  It  should  also  be  useful  in  assessing  the 
relative  merits  of  flexible  versus  rigid  linings. 

Since  the  method  of  analysis  is  intended  to  be  applicable  to 
structures  buried  in  soil,  an  inelastic  medium,  the  appropriate 
values  for  the  constants  relating  load  to  relative  displacement  and 
velocity,  the  virtual  mass,  and  etc.,  must  be  determined  experi¬ 
mentally.  A  test  program  to  investigate  the  problems  associated 
with  the  evaluation  of  the  parameters  for  cylindrical  structures 
under  static  load,  is  presented  in  Section  III. 
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2.  The  blast  induced  horizontal  displacements  of  the  free  field 
(  at  points  remote  from  ground  zero)  will  generally  be  nonlinear 
with  depth.  The  static  behavior  of  silo  linings,  as  a  result  of  the 
imposed  differential  horizontal  displacements,  is  investigated 
analytically  in  Section  IV.  Expressions  are  derived  for  the  determi¬ 
nation  of  bending  and  shear  stress  under  the  assumption  that 
loading  is  proportional  to  the  displacement  of  the  structure  relative 
to  the  given  free  field  profile.  It  is  shown  in  the  analysis  that 
axial  curvature  has  a  secondary  effect  on  the  shell  deformations  in 
the  plane  of  the  circular  cross  sections. 

3.  Theoretical  and  experimental  studies  concerning  the  axial 
skin  friction  loading  of  silo  linings  are  described  in  Section  V. 

These  surface  tractions  result  from  the  relative  vertical  displace  - 
ment  at  the  structure-soil  interface,  arising  from  soil  compres¬ 
sion  generated  by  the  intense  surface  overpressure.  The 
investigation  of  this  loading  aspect  relates  to  static  behavior  only. 
The  analysis  is  concerned  with  the  equilibrium  of  an  elemental 
wedge  of  soil  at  impending  failure  (bounded  by  assumed  failure 
surfaces  at  the  silo  lining  and  a  conical  surface  through  to  soil) 
assuming  Coulomb  shear  at  the  failure  surface.  Push-out  tests  of 
a  silo  model  embedded  in  dense  and  loose  Ottawa  sand  and  experi¬ 
ments  on  a  segmented  silo  model  in  dense  sand  are  also  given. 

4.  The  passage  of  the  air  shock  front  over  the  silo  cover  results 
in  the  sudden  application  of  vertical  load  to  the  silo  wall  if  the 
cover  is  integral  with  the  lining.  A  compression  wave  travels 
down  the  lining  and  is  reflected  at  the  base.  The  nature  of  the 
reflected  wave  depends  on  conditions  at  the  base.  In  general,  for 
silos  founded  in  soil,  the  elastic  modulus  of  the  foundation  will  be 
such  that  a  tension  wave  will  be  reflected.  Since  concrete  has  very 
little  tensile  strength,  the  vertical  reinforcing  steel  must  be  capable 
of  preventing  spallation  failures.  In  Appendix  A,  a  method  for  the 
analysis  of  the  effects  of  axial  stress  waves  is  presented  and 
applied  to  a  repr es entative  launch  silo. 
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THEORETICAL  ANALYSIS  OF  THE  RESPONSE  OF  SILO  AND 


TUNNEL  LININGS  TO  RADIAL  PRESSURE 

A.  Introduction 

The  stresses  and  displacements  produced  in  soils  by  the  air  shock 
of  a  nuclear  weapon  are  not  well  understood,  The  two  primary  reasons  for 
this  are,  first,  the  complexity  of  the  constitutive  equations  required  to 
describe  the  material  properties  of  soil  and,  secondly,  the  difficulty  in  the 
analysis  of  the  subsurface  effects  of  time-dependent  air-blast  loadings  in 
any  but  ideally  elastic  media.  The  semiempirical  methods  developed  for  the 
prediction  of  subsurface  stresses  and  displacements  as  functions  of  time  are 
based  on  greatly  simplified  analytical  models  and  the  limited  results  of  full 
scale  tests  (  see,  for  example,  Ref.  1)  . 

The  analysis  presented  in  this  section  gives  a  method  for  estimating 
the  significance  of  the  dynamic  effects  of  blast-induced  earth  pressure  load¬ 
ings  on  underground  silo  and  tunnel  linings  in  the  plane  of  the  circular  cross 
section.  It  is  not  an  exact  solution  (  in  the  sense  that  it  is  a  solution  for  the 
response  of  a  structure  in  an  inelastic  medium  v/ith  known  properties)  to  the 
disturbance  produced  in  such  a  medium  by  surface  air-blast  loading.  In 
fact,  no  exact  solutions  in  this  sense  have  been  obtained  for  the  design  of 
these  structures  for  static  earth  pressure  loadings.  That  is,  design  methods 
currently  used  for  static  loads  are  semiempirical.  It  is  hoped  that  with  the 
accumulation  of  experimental  data  on  the  response  of  buried  structures  to 
blast-induced  loadings,  the  analysis  presented  here  will  constitute  a  method 
of  design  for  dynamic  loads. 

B.  Basic  Assumptions 

The  objective  of  the  analysis  presented  in  this  section  is  the 
determination  of  the  significance  of  dynamic  effects  on  the  deformation  of  an 
originally  circular  cross  section  in  its  own  plane.  The  effects  of  variation 
in  the  loading  in  the  axial  direction  are  neglected.  This  is  equivalent  to  the 
assumption  that  the  loading  is  uniformly  distributed  along  the  cylinder  axis. 

In  this  way  the  problem  is  reduced  to  the  determination  of  the  response  of 
a  ring  to  forces  acting  in  the  plane  of  the  ring. 
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The  tangential  component  of  the  soil  pressures  acting  on  the 
ring  is  neglected.  The  magnitude  of  this  force  is  limited  by  the  coefficient  of 
friction  and  the  magnitude  of  the  normal  component  of  the  loading.  There 
does  not  seem  to  be  a  simple  way  to  estimate  these  forces. 

One  possible  approach  to  the  design  of  a  silo  or  tunnel  lining 
would  be  to  estimate  the  free-field  stresses  and  strains  using  methods 
currently  available  and  then  determine  which  of  the  following  conditions 
impose  the  most  severe  stresses;  1)  the  radial  component  of  the  free-field 
stress  acting  as  a  static  load  or  2)  the  stresses  caused  by  requiring  the 
radial  displacements  of  the  cylinder  to  conform  to  the  free-field  displace¬ 
ments.  These  two  conditions  represent  limiting  cases  for  static  loading.  In 

the  conventional  design  of  tunnels  for  static  loads,  this  procedure  has  been 

* 

refined  by  taking  into  account  the  flexibility  of  the  structure  (  3)  .  This  can 
be  done  by  assuming  that  differential  radial  displacements  of  points  on  the 
structure  relative  to  displacement  of  corresponding  points  in  the  free  field 
result  in  changes  in  the  loading,  directly  proportional  to  the  magnitude  of  the 
relative  displacement.  There  is  some  experimental  data  available  (3,  4)  , 
on  both  full  scale  and  model  structures  to  use  as  a  basis  for  estimating  the 
appropriate  value  of  the  constant  of  proportionality  (the  "foundation  modulus") 
for  various  types  of  soils. 

In  the  design  for  dynamic  loads  induced  by  air  blast,  it  is  clear 
that  effects  not  present  in  static  loading  may  be  important.  That  is,  one 
would  expect  the  loading  at  a  point  on  the  structure  to  be  influenced  by  the 
relative  velocity  as  well  as  the  relative  displacement  with  respect  to  a 
corresponding  point  in  the  undisturbed  free  field.  In  addition,  one  would  sup¬ 
pose  that  the  presence  of  the  surrounding  soil  has  the  effect  of  increasing  the 
natural  periods  of  vibration  of  the  structure.  An  effective  mass  has  been 
added  to  the  actual  mass  to  account  for  this  effect.  Since  there  is,  at  this 
time,  no  experimental  data  on  which  to  base  the  selection  of  appropriate 
values  for  the  constant  of  proportionality  between  relative  velocity  and 
change  in  load  or  the  virtual  or  effective  mass,  estimates  must  be  made. 

In  this  analysis,  an  attempt  is  made  to  take  each  of  these  effects 

$ 

Numbers  in  parenthesses  cite  references  listed  on  page  140. 

ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


-  6  - 


into  account.  We  assume  that  the  loading  on  the  structure  consists  of  the 
components : 

1)  the  radial  component  of  the  free-field  stress, 

2)  a  component  proportional  to  the  relative  displacement 

(  wq  -  w)  ,  where  w  is  the  radial  displacement  of  a  point  on  the 

shell  and  w  is  the  corresponding  component  of  the  free-field 

displacement,  and 

3)  a  component  proportional  to  the  relative  velocity  (w  -  w)  , 

where  dots  denote  differentiation  with  respect  to  time. 

A  virtual  mass  has  been  added  to  the  actual  mass  in  computing 
the  response  of  each  of  the  modes.  To  compute  the  loading  on  a  cylinder 
according  to  this  scheme,  it  is  necessary  to  know  the  free-field  stress  and 
displacement  as  functions  of  time  and  to  select  the  factors  of  proportionality 
which  relate  relative  displacement  and  velocity  to  load.  To  obtain  a  reasonable 
estimate  of  the  range  of  values  for  these  parameters,  the  properties  of  elastic 
media  (not  a  true  soil)  are  considered  (Section  II  E)  .  Ultimately,  the 
selection  of  these  factors  must  be  based  on  experimental  results  in  soil. 

As  mentioned,  methods  of  predicting  free-field  stress  and 
displacement  have  not  been  advanced  to  the  point  where  it  is  possible  to 
predict  these  quantities  as  functions  of  time  for  a  specified  range  of  blast 
and  soil  parameters.  It  is,  however,  possible  to  make  estimates  of  peak 
stress,  displacement,  and  strain.  To  estimate  the  importance  of  dynamic 
effects,  the  peak  stress  and  strain  in  the  plane  of  the  cross  section  are 
assumed  produced  by  the  passage  of  a  plane  compression  wave  traveling  at 
a  constant  speed.  The  results  should  permit  an  estimation  of  the  importance 
of  dynamic  effects. 

C .  Loading 

The  maximum  values  of  the  free-field  stress  and  strain  in  the 
plane  of  the  cross  section  as  well  as  the  rise  time  at  these  values  are 
assumed  known.  The  speed  with  which  the  disturbance  transits  the  cross 
section  is  assumed  known  and  constant.  For’  a  silo,  this  speed  will  be 
equal  to  the  speed  of  the  air-shock  front  since,  for  the  pressure  range  of 
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of  interest  (  above  100  psi)  ,  the  speed  of  the  air  shock  will  exceed  that  of 
compression  waves  in  most  soils.  For  a  tunnel  oriented  side-on  to  ground 
zero,  the  speed  with  which  the  cross  section  is  engulfed  cannot  exceed  the 
speed  of  compression  waves  in  the  soil,  while,  for  end-on  orientation,  less 
time  will  be  required  for  engulfment.  Sketches  of  each  of  these  cases  are 
shown  infigures  1  and  2.  In  the  case  of  the  silo  and  the  case  of  end-on 
orientation  of  the  tunnel,  the  speed  with  which  a  given  cross  section  is  engulfed 
is  greater  than  the  speed  of  a  compression  wave  in  the  soil.  Therefore,  the 
hypothetical  compression  wave  which  is  assumed  to  cause  the  observed  stress 
and  strain  in  the  plane  of  the  cross  section  must  make  an  angle  X  with  the 
normal  to  the  plane  of  the  cross  section  as  shown  in  the  sketches.  The 
assumption  that  the  observed  free-field  values  of  stress  and  strain  could  be 
produced  by  a  compression  wave  imposes  certain  restrictions.  Since  in  a 
plane  compression  wave  the  strains  in  the  plane  of  the  wave  front  are  zero, 
there  is  no  strain  in  the  direction  of  the  diameter  perpendicular  to  the  direction 
of  engulfment.  In  the  three  cases  shown,  this  diameter  is  perpendicular  to  a 
line  through  ground  zero  and  therefore  it  is  expected  that  in  these  cases  the 
observed  strains  will  be  negligibly  small.  However,  for  tunnel  orientations 
intermediate  to  those  shown  in  figure  2,  this  strain  will  be  somewhat  larger, 
but  it  will  doubtless  be  much  smaller  than  the  strain  in  the  vertical  direction. 

A  second  restriction  imposed  by  our  assumption  is  that  a  certain 
relationship  exists  between  the  observed  values  of  the  stress  biaxiality  ratio 
(  )  ,  the  apparent  modulus  (  E*  )  ,  the  ratio  of  stress  ( cr)  to  strain  (e)  in 

the  plane  of  the  ring,  and  the  speed  of  engulfment  (  c)  .  With  reference  to 
figure  l,  the  speed  of  compression  waves  for  a  given  soil  is 


where 


and 


(1) 


Ec  is  the  modulus  of  the  soil  observed  in  confined  compression, 


p  is  the  density  of  the  soil. 

If  the  value  of  c  is  specified,  the  angle  X  is  then 
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c 

~*-~Air  Shock  Front 


(b)  Plan 


Fig.  1  SILO  STRESS  WAVE 


Ground  Shock  F ront 
Elevation 

(a)  End  -On  Orientation 


c  =  c  /  cos^ 
c 

Section 


(b)  Side-On  Orientation 


Fig.  2  TUNNEL  STRESS  WAVE 


1C  - 


or 


where 


X  =  cos 


c 

c 

c 


X  cos 


-1 


2_ 

r 


Now  for  a  specified  value  of  the  biaxiality  ratio  (  jo )  the  stress 
in  the  direction  of  propagation  of  the  compression  wave  (  gJ  )  is  related  to 
the  stress  ( o~)  (  see  Fig.  i)  by  the  formula 


0-'  -  [/3  +  (  1  -  /3  )  r2  j  a-  , 

while  the  strain  («  )  is  related  to  the  strain  in  the  direction  of  the  wave 
(t  ')  by 

«  =  e  '/  r2 

The  apparent  modulus  (E  )  is  then 


E 


* 


G~ 


where  the  fact  that 


is  used. 


E 

c 


/ 3  +  (  1  -ft)  rl 


(V 


Equation  2  shows  that  it  is  not  possible  to  account  for  any  desired 
combination  of  observed  values  of  E^  ,  /z>,  and  c,  by  the  proposed  scheme. 
Given  any  two  of  these  quantities,  the  third  may  be  determined.  The  degree 
of  success  with  which  the  observed  values  can  be  fitted  by  the  proposed  method 
depends  upon  the  extent  to  which  they  actually  represent  the  effects  of  the 
passage  of  a  compression  wave. 

The  stress  wave  and  shell  coordinate  system  in  the  plane  of  the 
cross  section  are  shown  in  figure  3. 
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Fig.  3  STRESS  WAVE  AND  SHELL  COORDINATE  SYSTEM 

With  reference  to  the  coordinate  system  shown  in  the  sketch,  the 
variation  of  the  compressive  stress  in  the  x-direction  is  of  the  form: 

C~  =  q  ■  f(  x  -  ct)  for  x 

>  (3) 

0  ~  0  for  x  ^  ct  J 

The  compressive  stress  in  the  y-direction  is  assumed  equal  to 
0(3"  for  all  values  of  x,  that  is,  the  biaxiality  ratio  is  assumed  constant. 
The  strain  in  the  x-direction  will  be  assumed  proportional  to  the  stress. 
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The  displacement  in  the  positive 
originally  at  x  is  ,  for  small  strains, 


u(x) 


£  (x)  d  x 


x  -  direction  of  a  particle 


or  using  equations  3  and  4, 


u(x) 


f (x  -  c  t  )  d  x 


The  inward  radial  free-field  displacement  of  particles  originally  on  the 
surface  of  the  cylinder  is  then 

w  =  u  cos  9  ,  for  cos  9  ^  1  -  X 

o  ’ 


(5) 


where 


c  t 


Replacing  u  by  its  value  from  equation  5, 

'C  t 


w  =■ 
o 


C T  cos  9 
o 


*a(l  -  cos  ©) 


f  (x  -  c  t  )  d  x 


cos 


9  ^  1  -  X  (6a) 


w  *  0  ,  cos  9^-1  -X 
o  ’ 


(6b) 


The  inward  radial  velocity  of  a  narticle  is 

»c  t 


c O'  cos  9 
o 


■£ 


whe  re 


,  df 
f  (x  -  c  t  )  = 


d (x  -  c  t  ) 


f  •  (x  -  c  t  )  d  x 


cos  9) 


,  cos  0^1  -X 
(7a) 


w  =  0  ,  cos  9  ^  1  -  X 

o 


(7b) 
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The  radial  component  of  the  free- field  stress  is 


=  y—  ^  1  +/6  +  (1  -  fh  )  COS  2  ©J  ,  cos  0^1-t 
or  using  equation  3 

|”  l+y&  +  (l-/S)  cos  2  ej  f(x  -  c  )  ,  cos  ©  ^  1  -  "C 

(8a) 

rr  =0  ,  cos  ©  ^  1  -  T  (8b) 

urr 

Equations  6,  7,  and  8  give  the  radial  components  of  the  free- field 
displacement,  particle  velocity,  and  stress  in  terms  of  a  general  pulse  shape, 
f(x  -  c  t)  .  The  particular  pulse  shape  that  will  be  considered  in  this  study 
is  one  which  rises  linearly  to  a  constant  value  (see  Fig. 4  ).  This  pulse 
is  defined  as  follows! 

0,  x  -  ct 

f(x  -  c  t)  =  <  —  -  'J0-  -  y  c (t  -  t^)4x  4  c  t 
1  ,  X  c(t  -  t  ^) 


f(x) 


Fig.  4  PULSE  WITH  LINEAR  RISE 
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In  discussing  the  state  of  affairs  as  the  pulse  engulfs  the  cross 
section.it  is  convenient  to  define  the  following  angles  (See  Fig.  5) 


Fig.  5  WAVE  FRONT  POSITION  ANGLES 


f  -1/  , 

C08  (  1 

< 

TC 

V- 


r+  Tj)  .  rj*  -c  42+  tj 
.  2  +  Tj  4  -C 


fco  s  "* 1  (  l  -  t  )  ,  o  4  r  4  2 
[TC  ,  2  4  -c 


The  general  expressions  for  free-field  radial  displacement  particle  velocity 
and  stress  become, for  the  pulse  shape  being  considered, 

Case  1  (0  414^) 

for  -  0  4  0  4  0 


CTo 

w  *  - 

0  E* 

a  cos  0 

|  cos  0  -  (1 

(9) 

cos  0 

cos  0  -  (1 

->] 

(10) 

°  E* 

T1 
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(11) 


c r. 


rr 


0“o 

rr. 


£cos  0  -  (1  -  x)  J  £  1  +/?>  +  (1  -  fb  )  cos  20] 


b.  for  0  —  0  —  2  TC-  0 


w  ■=  w  =  cr  =  0 

o  0  rr 


Case  2  ("C^T  4  2  t-T^.) 


a. 


b. 


for  fi^e  4  ff  and  2  H  -  0  4  6  4  2TT  -  0X 

Equations  9,  10,  and  11  apply 
for  0^4©  4  2TT-  0 


w  «  w 
O  o 


0“  -  0 
rr 


for  -0  4  9  4  0 


(12) 


era 

r 

•Cl 

w  •=  ■ 
0 

0 

V* 

cos  0  cos  0  - 

u 

(1  -X)  -  — 

(13) 

E 

t 

• 

c  °"o 

w  «  ■ 
0 

E* 

cos  0 

(14) 

cr  - 

rr 

Co 

2 

[l  +  /2>  +  (1  - 

)  cos  2  ©j 

(15) 

Case  3  (  X  4  2  +  Xx) 


In  this  case,  equations  13,  14,  and  15  apply  for  all  values  of  0. 


The  Fourier  expansions  of  wq 


w  . 
0 


and  0"  have  the  form 
rr 


w 

o 


w 

o 


-I  ‘ 


cos  n  0 


cos  n  0 


(16) 


(17) 
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u  _ — 

^rr  "  ~T  +  2_  bn  cos  n  9 
n  -  1 


(18) 


where  the  coefficients  b  ,  c  ,  and  d  have  the  following  value  si 

n  n  n 

Case  1  (0 


b  - 


<T-a 


n  zrzrx 


<Toa 


n  e Vr 


-  (1  -  r  )  (1*0  )  s0  (0,  n)+  i-j^J  Sl  (0,  n) 

.  (1-g  )i<,x-Q  )  s?  (0(  n)+  ^  (0>  n),  (19) 

J-  (1  -  r  )  So  (0,  n)+  ±.  [--L+(l  -r)2]  Sl  (0,  n) 


-  (1  g—  S2  (0,  n)  +  -g—  S3  (0,  n) 


c 

d  -  - ° 

n  E *TTz- 


(20) 


SQ  (0,  n)  -  (1  -  T  )  Sx  (0,  n)  4-  S2  (0,  n)  l  (21) 


in  which 


S  (0,  n)  -  /  0 


S0  (0,  n)  -  0  ,  n  -  0 


and 


(22) 


S  <0,  n) .  ^  (-  - 1?  l. 4  sin  *  o  ?  ,  „  0 1 

l  n-i  n  +  i  * 


Si  (0,  „)  -  0  +  ^±jLL  ,  n-l 


i  -  1,  2,  3 
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Case  2  (  ^  4  2 +T1) 


bn  -  -^1  f  (1-yS)  So  (01>  n)  +  (0^  n)  -  (1  -  x)  (1+/3)  S'^  (n) 

S2'M  +  i*S  S3'  (n)  ]  ,  (23) 
So  (01  >  n)  -  [l  “  x  +  ]  S1  (01,  n)  +  -1-  S2  (01>  n) 


_  °~oa 
'n  TCe*  _ 

-  (1  -  X)  S  '(n)  + 


(24) 

^  +  (1  -  x)2  S1  (n)  -  — — S2  (n)-t--g-  S3  (i;)j 


n  TCe 


1S_|si  (0V  n)+  So'(n)  -  (1  -X)  ^'(n)*  -A-  S,,'  (n)j*  (25) 


where 


si-  ’  X 


l  L 


•Si  (0,  n)  -  Si  (01  ,  n) 


For  Xx  ^  0 


(n)  -  0  ,  For  X-L  =  0 
Case  3  2  +  T±) 


(26) 


3  -  cr  \ 

n  o  \ 


1  -*y2>  ,  n  -  0 
0  ,  n  ■  1 

(1  -  ft  )/2,  n  •=  2 
0  •»  n  -  3,  U, 


r 


c  = 
n 


^0  * 


1  ,  n  -  0 


^  (X  -  1  -  X1/2)  ,  n  -  1 
1/2  ,  n  =  2 

0  ,  n  =  2,  3,  . 


(27) 


(28) 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


18 


c  cr 


d  * 

n 


1 


0 

1  , 


,  n  -  0 
n  -  1 


(29) 


»  n  *  2,  3 > 


The  loading  normal  to  the  surface  of  the  cylinder  can  be  written 


as  follows j 


P  -  (Trr  +  k(wQ  -  w)  +  s(wq  -  w) 

in  which  k  and  s  are  factors  of  proportionality  relating  relative 
radial  displacement  and  velocity  to  load. 

Assuming  radial  displacement  of  the  cylinder  in  the  form 


(30) 


oo 


w  -  a 


n  *  1 


a  cos  n  0  , 
n  > 


the  load  acting  on  the  cylinder  can  be  expressed  as 
oo 


p  *  ^  p^  cos  n  Q 


n  *  0 


where 


and 


p  **  — i—  |  b  +  k(c  -  2  a  )  +  s  (d  -2a) 
^o  2l_°  0  0  0  ol 


p  =  b  ■+*  k  (c  -  a  )  +  s  (d  -  a  )  ,  n  **  1,  2, 
*n  n  n  n  n  n  ’  ’  ’ 


(31) 


"N 


J 


(32) 


The  coefficients,  ^>n  »  cn  >  and  equation  32  are  defined  in  equation  19 

through  equation  29  for  the  various  time  intervals. 
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D.  Equations  of  Motion 


The  radial  and  tangential  displacements  of  the  ring  are  assumed 


in  the  form 


w  •*  a  (t)  + 


n  =  1 


a  (t)  cos  n  0 
n 


z  ^ 


a  (t) 


sin  n  0 


n  =  1 


where  the  a  (t)  are  functions  of  time  which  are  to  be  determined.  The 
n 

equations  for  determining  these  functions  can  be  found  by  using  d'Alembert's 
principle  together  with  the  principle  of  virtual  work. 

The  intensity  of  the  inertia  forces  is,  in  the  radial  direction 


-  m  w  *  -  -7j—  a^  -  m 


r  00 

u  « 1  J 


cos  j  0 


and  in  the  tangential  direction 


-  m  v  «  -  m 


2Z  4-  3111  j®- 


j  - 1 


where  m  is  surface  density. 
Assuming  a  virtual  displacement. 


8 w  *=  8a  cos  n  0,  8 

n  ’ 


v  *  -  sin  n  0, 

n 


the  work  of  the  inertia  forces  is  found  to  be 


p2TT..  aS 

-  ra  a  /  w  Sw  d  0  -  m  a  / 

i)  J  ri 


v  6  v  d  0 


■=  -  2T[  a  m  a  8a 
o  o 


,  n  -  0 


-  IT a  mi£-+2L  V  Sa  ,n  *  0. 

y  n  n  ’  * 
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An  expression  for  the  strain  energy  of  a  thin  ring  of  thickness 


(h)  is 


Ta- 


(-3T  '  »>2  d  9 


ii  r.A 


( — +  w)  a  e 

d9 


where 


1  -  V 


Using  the  assumed  forms  for  w  and  v  given  in  equation  33, 


n.Iil,2  + 
a  0 


(nJ  -  l)2  a/  + 


n  *=  1 


9  _  oO 

Hj  ^  \  _  2  TCEI  \~  / 

u  — —  & +  — -?>  a0  +  —3—  2_(n 

12  a  2aJ  4 — 

n  -  1 


I  /  2  1  v  2  2 

"2_<n  '  1}  “n  • 


The  internal  work  done  on  a  virtual  displacement  is 
^  U  £ 

AW,  - x~— - -  6  a_  , 


therefore  , 


AW r  -  -  2TCEa  ±-  (1  +  -JL_ )  Sa  ,  n  -  0 
1  0  a  12  a2  0 
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A  ,,  TC  E  I  ,2  .  N  2  £ 

AW.  - - « -  (n  -  1)  a  o  a 

1  n  n 


n  f  0 


The  work  done  by  the  external  forces  is 
2ir 

SWQ«=a  /  p  S  w  d  ©, 


Sw  =  2TT  a  p  5a  ,  n  =  0 , 

O  r  A  A  '  * 


o  o 


and 


S  W  =  TC  a  p  S 


a  ,  n  4*  0. 
n  n  ' 


Now  from  the  condition  that  the  work  on  a  virtual  displacement  is  zero  it 
follows  that 

x2 

T 


-  2TC  a  m  a  -2"TCE  a  (1  +•  - )  2TC  a  p  =0 

o  o  a  ,  «  _  2  ro 


12  a 


or  replacing  p  by  its  value  from  equation  32  we  have 


m  a  +  s  a  +  j”  E  — 5  (1  +  - )  +  k 

L  a^  12  a 

for  the  equation  governing  pui-o  radial  motion. 
In  the  same  way  we  obtain 


a  *  — tt-  (b  +  kc+sd)  (37) 
o  2  o  o  o 


•  • 

2.m  a. 


ll  3  +  *  ax  c  *>]_  +  k  c1  +  s 

as  the  equation  governing  rigid  body  motion 
and 


m(n  +  1)  V  , 

— 5 — = - L—  a  +■  s  a  + 

n2  n  n 


E  I  ,  2  _.2  ^  , 

— r —  (n  -  1)  +  k 


(38) 


L  a 


a  =  b 
i  n  n 


+  k  c  +  s  d 
n  n 


(39) 
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as  the  equation  of  motion  for  n  ^  1  . 

In  the  last  of  the  above  equations  the  effect  of  circumferential  stress 
on  the  bending  deflection  has  been  neglected.  It  is  shown  in  5  that  in  the  case 
case  of  a  statically  loaded  ring  under  a  uniform  pressure,  pQ  ,  the  effect 
of  circumferential  stress  on  bending  can  be  taken  into  account  by  the  addition 
of  a  fictitious  pressure  equal  to 

oo 

-  ^f-  (^—!r  +  w)  ■  -f-  ^  (n2  -  !)  an  COS  n  e  .  (40) 

do  m 

n  -  2 

In  the  present  case,  dynamic  effects  will  be  ignored  in  determining  the  uniform 
component  (po)  of  the  radial  pressure.  This  simplification  uncouples 
the  equations  of  motion.  Then  from  equation  32, 

in  which  aQ  is  set  equal  to  its  value  under  static  loading, 

*0  *  "T  (V  k  “o'/  (1+  7775-  >  +  k]  (41> 

This  equation  is  obtained  from  equation  37  by  neglecting  the  inertia  and  damping 
terms.  Further  simplification  is  obtained  by  Introducing  the  approximation  that 


The  expression  for  pQ  then  becomes 

Po  "  ~T  (bo  +  k  °oV  (1  +  (42) 

E  h 

The  work  done  by  the  fictitious  pressure  on  a  virtual  displacement, 

Sw  *  S  a  cos  n  9 
n  ’ 

is 
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TCP.  (°2  -  !)  an  ^ 


an  "  * 


<bo  +  k  Co> 

2(1+ 

E  h 


(n2  -  1)  S 


n 


and  the  modified  equation  of  motion  becomes  for  n  ^  1 
2 


»<°7 12  v  -  3 » 


;  .  f  E 1 

a  +  — i - 

n  LIF^ 


E  1  (n2  -  l)2+k  an  -  bn 


+  k  c  4-  s  d 
n  n 


<V  k  °0) 

2(1+  ~2) 

E  h 


(n  -  1) 


n 


(43) 


Up  to  the  time  of  complete  engulfment,  that  is,  for  T  2  +  , 

the  equation  of  motion  for  the  bending  modes,  equation  43,  must  be  solved  by 

numerical  means.  After  complete  engulfment  the  b  ,  c  ,  and  d  are 

n  n  n 

constant  (equations  27,  28,  and  29),  so  that  solutions  in  closed  form  can  be  obtained. 

Xn  terms  of  the  dimensionless  displacement  and  time  variables, 

^  and  X  , 

x  *=  a  /a 
n  n 

X  -  ct/a 

Equation  43  becomes 


r2  xn"  +  fl  xn'  +  f0  xn  '  g-  "  '  2-  3- 


(44) 


where  primes  indicate  differentiation  with  respect  to  x  and  where 

2 


n2  +  1 


n 

c  s 


m  o 

Toa 


(45) 

(46) 
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(47) 


f 

o 


(n2  -  1) 

0“o 


(n2  -  1)  + 


k  a 


1 


O' 


2(1  + 


Eh  J 


and 


g  «  -j-  (1  ,  n  *  2 


g  -  0  ,  n  ^  2 


(48) 


The  character  of  the  solution  differs  depending  on  whether  or 
not  f Q  2k.  0  .  For  sufficiently  large  values  of  (^  ,  f  becomes  negative 
and  the  displacements  grow  without  bound,  buckling  occurs.  Since  yielding 
under  hoop  compression  occurs  before  the  buckling  load  is  reached,  this  case 
is  of  little  importance.  For  values  of  f  2k  0  the  solutions  are  either 
oscillatory  or  over  damped  depending  upon  whether  the  quantity  f^2  /U  f 
is  less  than  unity.  There  are  three  casesi 


Case  1  (f.  It  f  f_)  Under  damped 
_ X _ Q  <- _ 


k  -  J  (x  -  x  )  cos  ¥(  X  -2 -r,)-X-  f  x'  - 

n  y  n  ns  l  n 


•  v] 


sin  Y(z-Z-z. 


Case  2  (f  ^  h  f  f0)  Overdamped 
_ 1  o  ^ 


-2-Tl>x 


xn  ’  -b,  [(S*1))  K  -  x„s>  •  *’„] 

--K[<W  <xn  "  ’‘h  >  ‘  x,n]cC5*1l)<r"2'ri> 


Case  3  (f  t  “  U  f  f0)  Critically  Damped 

X  rs  c. 


(49) 


(50) 


V  ■  f(xn  -  xn  >  +Lx'n  ‘  $(xn  -  *„  >  ]  <  Z-Z-  T,>]  \ 

U  S  S  J  £ 


(51) 
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where 


s 


and 

x  «=  x  at  T  =  2  +  X, 

n  n  1 

x'  *  x'  at  »  2  +  X, 
n  n  1 

xu 

The  value  x  is  the  peak  static  displacement  of  the  n—  mode.  Since 
ns 

g  =  0  for  n  ^  2  the  first  bending  mode  (n  =  2)  has  the  only  non-zero 
static  displacement. 

The  critical  value  of  damping  at  which  the  motion  loses  its 
oscillatory  character  is  given  by  the  condition 
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The  term  involving  CT^  in  equation  52  derives  from  the  term  that  was  added  to 
equation  39  to  account  for  the  effect  of  direct  stress  on  bending.  When  this 
effect  is  negligible  the  value  of  scr  increases  with  increasing  mode  number. 
When  the  direct  stress  nears  the  buckling  stress,  the  quantities  in  the  brackets 
in  equation  52  reaches  a  minimum  for  a  value  of  n  2  and  scr  can  first 
decrease  then  increase  with  increasing  n.  In  any  case,  for  n  sufficiently 
large,  the  effect  of  damping  diminishes. 

E.  Selection  of  Parameter  Values 

Before  the  method  can  be  applied  to  the  investigation  of  the  dynamic 
effects  of  blast-induced  loads,  values  must  be  assigned  to  the  constant  k 
and  s  which  relate  changes  in  load  to  relative  displacement  and  velocity,  and 
values  must  be  assumed  for  the  virtual  mass. 

Static  load  tests  on  culverts  ranging  from  30  to  84  inches  in  dia¬ 
meter  indicate  that  the  foundation  modulus  (  k)  is  a  function  of  the  radius 
(  3)  .  Values  of  the  quantity  ka  where  a  is  the  radius  of  the  cylinder  seem  to 
be  constant  for  a  given  soil  condition.  Values  of  ka  ^  200  psi  are  cited  in 
reference  3  for  various  types  of  backfill.  Values  of  the  compression  moduli 
and  density  of  the  fill  material  are  not  given  so  it  is  difficult  to  relate  the 
observed  values  of  the  foundation  modulus  to  other  soil  properties. 

Watkins  (4)  obtained  values  of  ka  ranging  from  about  2300  to 
3700  psi  for  various  mixtures  of  silt  and  clay.  These  results  were  obtained 
from  tests  on  3-7/8-in.  diameter  cylinders.  Increasing  values  of  ka  were 
found  with  increasing  density  of  the  soil,  so  that  apparently  ka  increases 
with  increasing  values  of  the  compression  modulus. 

The  plane  strain  solution  for  the  displacements  produced  by  the 
application  of  a  radial  pressure  varying  as  cos  n  9  to  the  boundary  of  a  hole 
in  an  infinite  elastic  medium  is  geven  in  Appendix  B  as  follows 

n2  1 

ka  =  -%  "  t -  E  ,  n  k  2  (53) 

Thus,  for  an  elastic  medium,  the  foundation  modulus  increases  with  increas¬ 
ing  mode  number  and  is  directly  proportional  to  the  modulus  in  confined 
compression. 
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If  a  similar  solution  is  obtained  for  n  =  l  (rigid  body  displacement), 
the  foundation  modulus  is  found  to  be  zero.  That  is,  there  is  no  resistance 
to  the  displacement  of  a  hole  as  a  rigid  body  in  an  infinite  elastic  medium. 
For  this  reason,  a  zero  value  for  the  foundation  modulus  is  used  in  the 
computation  of  the  rigid  body  motion  (  n  =  1)  .  A  fixed  value  of  ka,  equal  to 
some  fraction  of  Ec,  is  used  in  the  computations  described  in  the  following 
subsection  for  all  other  values  of  n  since  there  is  no  basis  for  estimating 
what  the  variation  with  n  should  be  for  soils. 


There  is  no  experimental  data  on  which  to  base  the  selection  of 
the  value  of  s,  the  constant  of  proportionality  relating  changes  in  load  to 
relative  velocity.  The  choice  of  s,  determines  the  magnitude  of  the  jump  in 
pressure  at  the  head  on  point  (6=0)  at  t  =0.  There  is,  as  yet,  no  evidence 
that  such  reflection  effects  exist  in  soils;  however,  for  an  elastic  medium, 
the  reflected  pressure  would  be  twice  the  incident  pressure  for  the  case  of 
normal  incidence  of  a  compression  wave.  To  obtain  this  result  form  the 
an.\ lysis  presented  here,  a  value  of 


s  =  p  cc  (  Elastic  Medium) 


(54) 


would  be  required.  Selected  values  of  the  parameter,  s,  used  in  this  analysis 

are  taken  as  less  than  the  value  for  an  elastic  medium,  i.  e.  ,  - — <£.1.0. 

Pc 

c 

It  is  shown  in  reference  6  that,  in  the  analysis  of  the  response  of 
a  cylindrical  shell  submerged  in  an  incompressible  fluid,  the  interaction 
between  the  fluid  and  the  shell  can  be  completely  accounted  for  by  the  addition 
of  a  virtual  mass  to  the  actual  mass  of  the  shell.  This  virtual  mass  is  given 
by  the  formula 

n  -  (55) 


m 


=  P 


n 


^77 


where  p  is  the  density  of  the  soil,  a  is  the  radius  of  the  cylinder,  and  n  is  the 
displacement  mode.  It  seems  plausible  to  suppose  that  the  use  of  some 
fraction  of  this  virtual  mass  is  appropriate  for  an  underground  structure.  In 
the  computations  described  in  the  following  section,  a  mass  equal  to  K  times 

y 

m  is  added  to  the  actual  mass  of  the  shell,  where  the  factor  K  is  between 
v  '  m 

n  ,  v 

zero  and  unity. 
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F.  Results  of  the  Analysis 

It  is  clear  that  the  utility  of  the  proposed  method  for  the  analysis  of 
the  response  of  silo  and  tunnel  linings  to  blast-induced  loadings  can  be  established 
only  by  demonstrating  agreement  between  predictions  of  the  method  and 
experimental  results.  There  is,  at  this  time,  no  data  that  can  be  used  for 
this  purpose.  It  is  expected  that  tests  will  be  carried  out  in  the  near  future 
to  check  the  theory. 

In  the  interim,  something  is  to  be  gained  by  an  investigation  of  the 
nature  of  the  solutions  given  by  the  method  and  by  the  study  of  the  effects  on 
these  solutions  of  variations  in  the  parameters  whose  values  at  this  time  are 
more  or  less  indeterminate.  Although  it  is  not  possible  to  assign  precise 
values  to  the  parameters  on  which  the  solution  depends,  it  may  prove  that  for 
wide  variations  in  their  values  the  stresses  remain  within  certain  limits.  In 
this  way  it  may  be  possible  to  draw  useful  conclusions  concerning  the 
significance  of  dynamic  effects.  In  other  words,  it  may  be  possible  to  show 
that  with  the  application  of  an  appropriate  factor  to  account  for  dynamic  effects, 
an  analysis  for  static  loading  is  adequate  for  design. 

A  further  reason  for  investigating  solutions  given  by  the  method  is 
that,  in  this  way,  a  fairly  accurate  knowledge  of  the  relative  effects  of  vari¬ 
ations  in  the  flexibility  of  the  structure  can  be  obtained.  The  effects  of  vari¬ 
ation  in  wall  thickness  in  a  concrete  cylinder  and  variations  in  the  parameter 
(s)  were  considered  in  the  first  interim  report  on  this  program  (  7)  .  However, 
the  results  suffer  from  the  fact  that  the  value  of  the  apparent  modulus  E  used 
in  the  three  cases  considered  did  not  satisfy  the  requirement  contained  in 
equation  2.  Since  too  low  a  value  of  E  was  used,  the  results  exaggerate  the 
dynamic  effects. 

In  the  present  report,  typical  examples  of  curves  of  displacement, 
pressure,  and  stress -vers us -time  are  presented.  Because  these  are  the 
quantities  that  will  be  measured  in  any  future  experimental  investigations,  it 
is  of  interest  to  note  the  characteristics  of  the  response  predicted  by  the 
theory.  No  attempt  has  been  made  to  study  completely  the  effects  of  variations 
in  the  parameters  entering  into  the  solution. 
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Solutions  were  obtained  for  31  cases  listed  in  table  1  which  are 
intended  to  be  representative  of  a  tunnel  oriented  side-on  to  the  direction  of 
the  compression  wave  propagation  (  see  Fig.  2)  .  Solutions  for  seven  additional 
cases  listed  in  table  2  that  are  representative  of  a  silo  (  see  Fig.  1)  have 
also  been  obtained.  Six  displacement  modes  (n  =  0  through  5)  have  been 
included  in  the  solutions.  The  numerical  integration  was  performed  on  the 
Foundation's  UNIVAC  1105,  Remington  Rand  Computer.  The  computer 
program  is  given  in  Appendix  C. 

1.  Tunnel  Cases  (  Side-On  Orientation) 

Variations  in  the  values  of  ka,  the  virtual  mass,  s,  Z 7^  and  the 
ratio  of  radius  to  thickness,  a/4i  were  considered.  The  following  fixed  values 
were  used  for  the  other  constants. 


200  psi 
l/3 

=  21,  600  psi 

cc  =  1 , 000  fps 

-4  2-4 

1.  5  x  10  lb-sec  -in.  ,  soil  density  (  100  pcf) 

-4  2  -4 

2.25  x  10  lb-sec  -in.  ,  density  of  shell  material 

(  concrete) 

3  x  10^  psi  (  appropriate  for  concrete) 


The  compression  modulus  is  in  the  range  of  values  for  representative 
soils  and  a  biaxiality  ratio  (  )  equal  to  one -third  is  justifiable  for  the  side- 

on  orientation  since  this  value  of  lateral  pressure  to  applied  normal  pressure 
is  commonly  used  for  a  wide  variety  of  soils. 


Three  values  of  ka  equal  to  5000,  2000  and  500  psi  were  used. 

These  variations  cover  the  extremes  of  the  range  of  values  given  in  references 

3  and  4.  Values  of  the  quantity - equal  to  0,  0.  125,  0.  25,  and  0.  5  wfere 

P^c 

used  to  investigate  the  effects  of  variations  in  s  while  values  of  Km  equal  to 
zero  and  one-half  were  used  to  study  the  effects  of  variations  of  virtual  mass. 
Values  of  a/h  equal  to  30,  20  and  10  were  used  to  examine  the  effect  of  thick- 
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ness  variations.  These  ratios  represent  wall  thicknesses  of  6,9  and  18 

inches  respectively  for  a  tunnel  radius  of  15  feet.  Constant  values  of 

s/oc  equal  to  0.  125  and  K  of  0.  5  were  selected  for  the  a/h  variation. 

'  r  c  m 

The  majority  of  the  solutions  were  obtained  for  a  step-pulse  loading  (77^  =  0)  ; 
however,  the  effect  of  a  rise  time  equal  to  one  transient  time  of  the  stress 
wave  {~Cy  =  2)  was  considered  for  each  ka  with  fixed  values  of  =  l/ 2, 
s  /pcc=  l/8  and,  a/h  =  10.  For  a  tunnel  radius  of  1  5  feet  the  rise  time,  t^, 
equals  30  msec. 

Typical  curves  of  pressure,  relative  displacement,  and  maximum 
combined  stress-versus-time  are  given  in  figure  6a,  b  and  c.  The  peak 
values  of  these  quantities  and  the  times  of  their  occurrence  are  listed  in 
table  1.  Curves  of  velocity-versus-time  for  rigid  body  response  are  shown 
in  figure  7. 

2.  Silo  Cases 

Two  values  each  of  ka,  Q-  and  c  were  considered.  A  concrete  lining 
with  ifcsy.o  ratios  of  radius  to  thickness  and  a  single  steel  lining  were  studied. 
Fixed  values  of  constants  for  all  cases  were 

K  =  0.  5  c  =  1000  fps 

m  c  r 

v 

Tj  =  0.0  s/pc^  =  0.25 

p  =  1.5x10^  lb  -  sec^  in  ^  /3  =1.0 


and 


Ec  =  21, ,600  psi 


^s 


E 


I"  2.  25  x  10  ^  lb  -  sec^  -  in  ^ , 
[_  7.  35  x  10  ^  lb  -  sec^  -  in 
C  6 

3  x  10  psi,  concrete  shell 
^30  x  10^  psi,  steel  shell 


concrete  shell 
steel  shell 
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Fig.  6c 


Nondimensional  Rigid  Body  Velocity, 


Fig.  7  VELOCITY -TIME  VARIATION  (TUNNEL) 


Fig.  8a  PRESSURE -TIME  VARIATION  (CASE  35) 
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Fig.  8c  COMBINED 


I 


nsional  Time,  ~C 


5-TIME  VARIATION  (CASE  35) 


Nondimensional  Rigid  Body  Velocity 


Nondimensional  Time,  T7 
Fig.  9  VELOCITY-TIME  VARIATION  (SILO) 
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A  biaxiality  ratio  of  unity  has  been  chosen  since  the  plane  of  the 
ring  for  a  silo  lining  is  horizontal  and  the  free-field  pressures  in  the  x  and  y 
directions  (  see  Fig.  1)  are  both  normal  to  the  air-blast  overpressure  loading 
acting  on  the  ground  surface.  The  apparent  soil  modulus,  E*,  has  been  based 
on  equation  2  and  is  determined  once  the  engulfment  speed,  c,  is  selected. 

For  an  air-blast  overpressure  of  300  psi,  a  value  of  c  =  4810  fps  is  obtained 
for  a  standard  sea  level  atmosphere.  This  value  is  used  for  the  four  cases 
(  cases  32  through  35)  selected  to  study  the  effect  of  CTq  and  ka  variations. 

An  engulfment  speed  (  c)  of  2000  fps  was  used  to  investigate  the  effects  of 
lining  thickness  variations  (Cases  36,  37  and  38)  to  accentuate  the  dynamic 
effects.  This  engulfment  speed  grossly  underestimates  the  speed  associated 
with  the  peak  free-field  stress  of  200  psi  considered  in  the  thickness  variation 
study  since  an  air  blast  overpressure  of  only  37.  5  psi  produces  an  air  shock 
front  velocity  of  2000  fps. 

Typical  curves  of  pressure,  relative  displacement,  and  combined 
stress-versus-time  are  given  in  figures  8a,  b  and  c  and  rigid  body  velocity- 
versus-time  curves  in  figure  9.  The  peak  values  of  the  response  quantities 
and  the  time  of  occurence  are  listed  in  table  2. 

G.  Discussion  of  Results 

The  ratios  of  maximum  dynamic-to- static  values  (  D.  F.  L.. )  of 
peak  combined  stress  and  relative  displacement  are  given  in  tables  1  and  2  for 
the  tunnel  and  silo  cases  respectively.  The  tunnel  cases  show  that  the 
introduction  of  a  relative  velocity  component  greatly  reduces  the  maximum 
response.  In  the  cases  where  ka  =  500  psi,  the  D.  L.  F.  for  stress  is  reduced 
from  approximately  1.8  for  zero  damping  to  a  point  where  the  static  values  are 
not  even  attained  during  five  transit  times  (  Z  —  10)  when  the  value  of 
s/pcc=  0.  5.  It  is  also  noted  that  even  relatively  small  values  of  s / pc  result 
in  large  D.  L.  F.  reductions.  The  effect  obtained  by  increasing  s/pcc  for  all 
cases  studied  corresponds  to  the  belief  that  the  response  of  buried  structures 
is  heavily  damped.  Since  all  solutions  are  for  constant  peak  free-field  stress 
of  O'  ,  the  effects  of  a  decaying  pulse  will  be  considerably  less  if  the  decay 
is  substantial  in  four  or  five  transient  times. 

The  value  of  s/pcc  determines  the  reflected  pressure  at  the  head-on 
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point  at  Z  —  0,  and  reflection  factors  accompanying  the  heavily  damped 
response  are  small  enough  to  seem  realistic.  When  the  pulse  has  a  rise  time 
of  =2.0,  the  reflection  effects  are  nearly  eliminated. 

The  effect  of  adding  virtual  mass  is  to  increase  the  maximum  response 
for  the  step  pulse  considered.  The  increase  in  the  D.  L..  F.  is  less  than  0.26 
for  all  cases  when  s/pcc^0.  125.  Mass  effects  would  be  much  more 
pronounced  for  more  realistic  loading  (decaying  free-field  stress  wave). 
Decreasing  the  wall  thickness  (  increasing  a/h  from  10  to  20  then  30) 
decreases  the  dynamic  effects  for  all  values  of  ka  considered.  The  reduction 
in  D.  L.  F.  is  less  than  0.  28  for  values  of  ka  equal  to  5000  and  2000  psi; 
however,  a  very  pronounced  effect  is  noted  when  ka  =500  psi.  It  is  observed 
that  the  static  relative  displacement  and  stress  enlarge  to  intolerable  values 
for  both  a/h  increases;  however,  the  static  values  are  not  obtained  within 
five  transit  times.  The  large  response  times  associated  with  these  cases 
could  conveivably  result  in  acceptable  displacements  and  stress  for  a 
decaying  free-field  stress  pulse. 

For  the  silo  cases  considered,  the  peak  free-field  stress  variations 
(Cases  32-35,  Table  2)  slightly  increase  the  dynamic-to-static  values  of 
maximum  relative  displacement  and  have  no  effect  on  the  combined  stress 
D.  L.  F.  All  absolute  dynamic  displacements  are  low  for  concrete  linings, 
even  for  the  large  D.  L.  F.  's  associated  with  Cases  34-36,  since  the  static 
displacement  is  very  small  for  reasonable  lining  diameters.  Bending 
stresses  are  almost  trivial  for  the  concrete  lining  cases.  These  low  values 
can  be  attributed  to  the  hydrostatic  free-field  stress  ( =  1)  and  the  large 
apparent  soil  modulus  and  engulfment  speed.  Even  the  unrealistic  reductions 
of  c  and  E*  considered  in  Cases  36  and  37  do  not  generate  appreciable 
displacement  or  bending  stress  response  values.  When  a  steel  shell  with  a/h 
of  60  is  assumed  (Case  38)  static  solutions  are  intolerable  (  e.  g.  ,  static 
combined  stress  of  775.9  x  200  —  155  ksi)  ;  however,  the  dynamic  relative 
displacement  and  bending  stress  are  less  than  0.  01  times  the  static  values 
up  to  10  transit  times  (t  =  20)  .  Again,  just  as  it  was  noted  for  the  tunnel 
cases,  the  large  response  times  associated  with  flexible  linings  will  have  an 
appreciable  effect  on  the  response  values  for  a  decaying  pulse. 
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The  curves  of  velocity  in  the  rigid  body  mode,  given  in  figure  7 
for  the  tunnel  study,  indicate  that  in  cases  where  the  virtual  mass  is  zero, 
the  velocity  overshoots  the  free-field  value.  The  addition  of  virtual  mass 
results  in  what  appears  to  be  more  "realistic"  behavior,  that  is,  the  velocity 
of  the  structure  approaches  the  free-field  value  asymptotically.  However, 
it  remains  for  experimentation  to  show  what  "realistic"  behavior  is. 


H.  Solution  for  Static  Loading 


The  results  of  the  preceding  section  show  that  in  the  cases  considered 
the  maximum  dynamic  stress  exceeds  the  stress  under  static  load  by  less 
than  40  percent  when  damping  of  s/pcc  ^  0.  125  is  included.  In  the  dynamic 
case  a  numerical  integration  of  the  equations  of  motion  is  required  for  the 
determination  of  the  displacement  amplitudes,  x^;  however,  for  the  static 
case  they  are  given  by  the  simple  formula 


Xn  =  g/fo’  n  =  2’  3> . > 

where  equations  47  and  48  define  f  and  g  respectively.  Since  g  =  0  for 
n  is.  2  there  is  no  bending  in  the  higher  modes  (n  ^  2)  .  The  maximum 
combined  stress  is  therefore  equal  to  the  sum  of  the  direct  stress,  f^,  plus 
the  bending  stress  in  the  first  bending  mode  (  n  =  2)  ,  f^. 


The  direct  stress  is  equal  to 

cos  n  9,  n  =•  0, 

1 

The  value  pQ  is  given  in  equation  42.  From  equation  32, 
effects, 


2,  3.  .  .  (  56) 

ignoring  dynamic 


p  =  b  4-  ka(c  -  x  )  ,  n  ^  1 
rn  n  v  n  n'  ’ 


Substituting  from  equations  27  and  28  and  noting  that  x^  = 


-(  14/34- ka/E*) 


i 

i 


0  for  n  2  we  find 
.  ka 

+  “ZE*  + 


ka 


°"o 


cos  29, 


(57) 
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where  positive  direct  stress  denotes  tension. 

The  displacement  ratio  in  the  first  bending  mode  is  obtained  by 
neglecting  the  dynamic  terms  in  equation  44,  thus 


xz=e/f0 

Using  equations  47  and  48  this  may  be  expressed  as 


or 


x2  = 


Go 

~r 


1  -  /S  +  k  a/E* 


E  I  k  a 

3  — -  4-  —m —  -  cr 


1  +  fb  +  ka  /  E* 
°  2(i  +  ka  V  "Eh) 


°S  (>+|-)  (iyJ  +  ^)/(1+/S+^  ) 


x2= 


2(1+  ka2/  Eh) 
l  +  y3  +■  ka/fe* 


,  El  ,  ka 

3  — 3-  +  3- 


°o 


(58) 


The  denominator  in  equation  58  becomes  zero  when  reaches  the  value 
required  to  cause  buckling  in  the  mode  n  =  2 

2 


iSJL  \ 

n,  K  Eh  '  /  ,  E  I  ka, 

CT  _  =  - ^ -  (  3  — ij—  +  -j)  . 


cr. 


Equation  58  can  also  be  written 

cx. 


X2~  1  -  ot- 


1-/a  +  F* 

T777^7H|T 

E 


Eh 


(59a) 


(59b) 


where 


oc 


2 


(59c) 


The  buckling  pressure  0“  is  so  large  for  reasonable  values  of 

cr2 

ka  that  the  effect  of  direct  stress  on  bending  in  the  second  mode  can  generally 
be  neglected.  In  this  case  <.  <  1  and  the  approximation  that 
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(59d) 


fD  +  fB 

Oo  J  G  =  0 

for  four  values  of  the  biaxiality  ratio  (  ^  )  and  values  of  the  radius  to 
thickness  ratio  (  a/ h)  from  10  to  30.  The  second  order  contribution  of  the 
bending  mode  (  n  =  2)  to  the  direct  stress  has  been  neglected  in  the  calcula¬ 
tions  (i.e.  ,  the  second  term  in  the  right  hand  member  of  equation  57  is 
neglected)  .  Values  of  E  =  3  x  10  psi,  E*=ka  and  four  values  of  the  founda¬ 
tion  modulus  (  ka  =  500,  2000,  5000  and  20,000  psi)  have  been  used  in  the 
computation. 

The  solid  curves  of  figure  10  represent  the  nondimensional 

combined  stress  when  the  effects  of  axial  ring  force  on  bending  are  neglected 

(  i.  e.  ,  utilizing  the  approximation  for  f^given  in  Eq.  61)  .  These  curves 

give  good  approximations  to  the  combined  stress  when  the  buckling  pressure, 

0"  ,  is  large  compared  to  o~ ,  or  when,  o£?=  O'*  /<J~  —  0.  The 

cr,  o  c  o /  ct  £ 

dasned  curves  give  the  exact  combined  stress  for  an  applied  pressure  of 
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q~o  =  200  psi  and  indicate  the  error  involved  in  the  approximation.  Observe 
that,  for  large  values  of  the  foundation  modulus,  the  approximate  solution 
yields  nearly  the  same  combined  stress  as  the  exact  solution  for  the  para¬ 
meters  selected. 

The  effects  on  maximum  combined  stress  of  increasing  the  wall 
thickness  while  holding  the  radius  fixed  can  be  seen  from  the  curves  in 
figure  10.  These  results  show  that  maximum  combined  stress  can  actually 
be  increased  by  increasing  the  wall  thickness  (decreasing  a/h)  Further 
examination  of  the  curves  reveals  that  variations  of  the  radius  to  thickness 
ratio  from  10  to  30  have  a  relatively  minor  effect  on  the  combined  stress 
for  2000^1ka  ^.5000.  As  expected,  increasing  the  biaxiality  ratio  and  the 
foundation  modulus  give  reductions  in  combined  stress.  However,  when 
ka  =20,  000  psi  and  a/h  ^15,  increasing  the  biaxiality  ratio  actually 
increases  the  combined  stress.  This  occurs  because  the  bending  stress 
becomes  relatively  unimportant  for  sufficiently  high  values  of  ka  and  ft  and 
the  direct  stress  (which  increases  with  3)  dominates  the  combined  stress 
expression. 
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Static  Combined  Stress  Ratio 


Radius  to  Thickness  Ratio,  a/h 

Fig.  10a  STATIC  COMBINED  STRESS- RADIUS  TO  THICKNESS 
RATIO  VARIATION 
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SYMBOL  INDEX 


Section  II 


a 

shell  radius 

ll 

rise  time  of  stress  wave 

a 

n 

peak  shell  modal  displacement 

u 

free -field  displacement 

bi 

Fourier  series  coefficients  of  <rrr 

U 

strain  energy 

c 

engulfment  velocity  of  shell  section 

V 

shell  tangential  displacement 

cc 

soil  compression  wave  velocity 

w 

shell  radial  displacement 

c. 

1 

Fourier  series  coefficients  of  w 

o 

w 

o 

free -field  radial  displacement  at  shell  surface 

d. 

1 

Fourier  series  coefficients  of  w 

0 

We'Wi 

external  and  internal  work 

£ 

Young's  modulus  of  shell 

x,y 

coordinates 

E 

e/(i  -  V2) 

X 

n 

a  /  a 
n 

£ 

apparent  modulus  of  soil 

*2 

Ec 

soil  modulus  in  confined  compression 

free  -  field  stress  biaxial  ratio 

VVf2 

equation  of  motion  coefficients 

e 

free  -field  strain  in  x-direction 

fD 

shell  static  direct  stress 

£o 

peak  free -field  strain  in  x-direction 

*B 

shell  static  bending  stress 

€' 

free-field  strain  in  direction  of  stress  wave  propagation 

g 

equation  of  motion  coefficient 

closed  form  solution  constants 

h 

shell  thickness 

o 

angular  coordinate  of  point  on  shell 

k 

foundation  modulus 

\j 

Poisson's  ratio 

K 

m 

v 

virtual  mass  coefficient 

P 

soil  density 

m 

V 

n 

virtual  mass 

shell  density 

n 

shell  displacement  mode 

<r 

free-field  pressure  in  x-direction 

P 

pressure  normal  to  shell  surface 

^o 

peak  free-field  pressure  in  x-direction 

Pn 

peak  modal  normal  pressure 

rcr2 

shell  buckling  pressure  in  second  mode 

Po 

uniform  radial  pressure 

^rr 

radial  free-field  pressure  at  shell  surface 

r 

c/cc 

0-' 

free-field  pressure  in  direction  of  stress  wave  propagation 

8 

coefficient  of  damping 

X 

ct/a 

8cr 

critical  damping  coefficient 

ctj/a 

S.,  S'. 

V  1 

t 

Fourier  series  arguments 

time 

a,«i 

X 

angular  coordinates 

angle  between  shell  axis  and  stress  wave  front 
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III.  STATIC  LOAD  TESTS  ON  TUNNELS 


A.  Introduction 

As  part  of  the  investigation  of  silo  and  tunnel  linings,  a  series 
of  static  tests  on  6-inch  diameter  tubes  (horizontally  oriented)  buried  in 
dense,  dry  Ottawa  sand  were  performed.  The  objective  of  the  test  program 
was  to  investigate  the  ability  of  the  analysis  given  in  the  Section  II  H  (  Static 
Response  of  Silo  and  Tunnel  Linings  to  Radial  Pressure)  to  predict  stresses 
and  displacements  under  static  load.  These  experiments  are  a  necessary 
preliminary  to  dynamic  tests  and  also  serve  to  point  out  some  of  the 
problems  associated  with  small  tunnel  model  testing. 

Tests  were  performed  on  steel  cylinders  with  radius  to  thickness 
ratios  of  60  and  40  (  one  each)  and  an  aluminum  cylinder  with  a  ratio  of  60 
in  a  sand-filled  pressure  vessel  three  feet  in  diameter  and  height.  The  test 
setup,  test  results  and  discussion  of  the  results  with  respect  to  the  theory 
are  presented  in  the  following  sections. 

Tunnel  model  strains  and  radial  displacement  measurements  were 
taken  near  the  center  of  the  model.  Instrumentation  for  the  tests  included 
both  electronic  and  mechanical  systems.  Electronic  measurements  were 
made  of  radial  displacements,  and  bending  and  direct  strain  in  the  tunnel 
lining.  Mechanical  measurements  were  made  to  ascertain  applied  pressures, 
soil  surface  deflections  and  rigid  body  tunnel  displacements. 

B.  Test  Setup  and  Instrumentation 

All  tests  of  the  tunnel  models  were  performed  in  a  3-ft  inside 
diameter,  3-ft  high  flanged  steel  cylinder  (  Fig.  11)  .  The  inside  surface  of 
the  l/2  -in.  cylinder  wall  was  lined  with  two  unbonded  sheets  of  Teflon 
(0.01-in.  thickness)  to  reduce  vertical  soil  restraint  at  the  containment 
vessel  wall  to  a  tolerable  value.  Experiments  (  8)  indicate  that  such  a  treat¬ 
ment  is  worthwhile  even  though  it  does  not  completely  eliminate  wall  friction. 
Cover  plates  2-1/8  in.  thick  bolted  to  the  cylinder  flanges  with  thirty-six 
1/2  in. -diameter  high-strength  bolts  and  a  rubber  sealing  membrane  placed 
over  the  sand  surface  and  between  the  cylinder  flange  and  upper  cover  plate 
complete  the  containment  vessel.  Surface  loading  is  accomplished  by 
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,  Sealing  Membrane 


Strain  Gage 


Section  1 


Fig.  11  TUNNEL  MODEL  TEST  ARRANGEMENT 
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admitting  air  pressure  between  the  membrane  ar.d  the  cover.  Two  diametrically 
opposite  holes  (  7-in.  diameter)  are  located  at  mid-depth  in  the  cylinder  wall. 
The  tunnel  model  extends  through  these  openings  and  projects  approximately 
l-i/2  in.  past  the  wall  at  each  opening.  Leather  boot  seals  were  provided  at 
the  model-containment  vessel  interface  to  prevent  soil  loss  in  a  manner  that 
did  not  appreciably  restrict  the  model  or  soil  motions.  A  sketch  of  the 
pressure  vessel  with  tunnel  model  in  place  3s  shov/n  in  figure  11.  An  overall 
view  of  the  test  setup  is  shown  in  figure  12a  and  the  model  just  prior  to 
burial  in  figure  i2b. 

All  models  were  cold-finished  from  cold  rolled-sheet  stock.  Two 
models  were  fabricated  from  mild  steel  with  radius  to  thickness  ratios  of 
40  and  60  (  0.  075  and  0.  05  in.  wall  thicknesses  respectively)  .  The  third 
model  was  6061-6T  aluminum  with  a  yield  stress  cf  35,  000  psi  and  radius 
to  thickness  ratio  of  60  (  0.05  in.  wall  thickness) 

All  strain  gages  used  for  the  test  series  were  temperature-com¬ 
pensated  and  epoxy  backed.  A  sing'e  layer  of  Scotch  tape  was  placed  over 
each  outer  gage  for  protection  purposes.  Bending  and  direct  stress  were 
determined  from  strain  measurements  made  witl  Budd  etched  foil  strain 
gages,  type  C6-141.  These  gages  are  produced  from  thin  alloy  foil  sheet 
by  photo-etching  techniques  similar  to  those  used  in  producing  printed  circuits. 
The  gage  is  designed  to  be  sensitive  to  the  rormai  stress  in  steel  in  the 
direction  of  the  gage.  The  grid  length,  width  and  the  overall  gage  width  are 
0.  25  in.  and  the  overall  gage  length  ’S  0.  503  in.  Gages  were  located  in  pairs 
on  the  inner  and  outer  surface  of  the  shell  to  measure  hoop  strain.  For  each 
of  the  three  tunnel  lining  models  gage  pairs  were  located  in  the  invert,  one 
spring  line,  45  degrees  above  the  spring  line  ar.d  the  crown  as  shown  in 
Section  1,  in  figure  il.  Strain  gages  are  also  mounted  at  Section  2  on  the 
outside  surface  at  the  crown  and  invert  to  measure  longitudinal  strain.  A 
J3aldwin  SR-4  battery  operated  strain  indicator  was  used  to  obtain  readings 
of  each  gage  for  all  surface  pressures. 

The  radial  displacements  of  the  tunnel  model  were  measured  at 
the  Section  2  location  shown  in  figure  li.  A  sketch  of  the  device  used  for 
these  measurements  is  shown  in  figure  13.  The  displacement  gage  consists 

of  a  Schaevitz  differential  transformer  with  a  l/2--in.  range  mounted  normal 
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Fig.  12a  TUNNEL  MODEL  TEST  SETUP 


12b  TUNNEL  MODEL  PRIOR  TO  SAND  PLACEMENT 


to  the  axis  of  a  2-in.  outside  diameter  (  10  B.  W.  gage)  steel  pipe  supported 
at  the  ends  of  the  tunnel  lining  model  by  spherical  bearings.  A  pointer  was 
attached  to  one  end  of  the  pipe  to  indicate  the  position  of  the  displacement 
probe.  Adjusting  screws  were  provided  at  the  bearing  assembly  to  initially 
center  the  rotation  axis  at  the  displacement  measurement  station.  Figure  14 
shows  three  views  of  the  radial  displacement  gage  and  the  Daytronic 
amplifier- indicator  used  to  indicate  the  displacement  reading. 

Displacements  of  the  upper  soil  surface  were  obtained  for  each 
test  by  dial  gage  readings  that  indicated  the  displacement  of  four  threaded 
rods  mounted  to  the  pressure  vessel  cover  as  shown  in  figure  15.  The 
threaded  rods  were  screwed  down  until  contact  of  the  rod  with  a  small  brass 
disk  cemented  to  the  sealing  membrane  completed  a  simple  electrical  light 
circuit.  The  rods  were  located  on  a  ray,  normal  to  the  tunnel  axis,  emanating 
from  the  center  of  the  pressure  vessel  cover.  Gages  were  positioned,  zero, 

2,  6,  and  10  in.  from  the  center.  Dial  gages  were  also  used  to  obtain  the 
rigid  body  tunnel  displacement.  A  gage  was  attached  to  the  containment  vessel 
at  the  points  where  the  tunnel  model  protruded  and  measurements  were  made 
of  the  radial  displacement  gage  support  plate  (which  was  rigidly  attached  to 
the  end  of  the  tunnel  model),  (  Fig.  14b)  . 

The  soil  used  for  all  tests  was  dense  Ottawa  sand  with  grain  size 
distribution  as  in  figure  16.  A  one-horsepower,  iO,  000-rpm  concrete  vibrator 
was  used  to  obtain  an  average  density  of  109.  0  pcf  for  each  test. 


The  soil  placement  procedure  was  as  follows: 


Step  1 

Soil  was  placed  to  the  tunnel  invert  (  15-in.  lift)  and  compacted  by 
penetration  of  the  vibrator  probe  to  the  lift  depth  in  the  following  pattern: 


Tunnel  Model 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOIOGY 


-  58  - 


Fig.  14a  RADIAL.  DISPLACEMENT  Fig.  14b  DISPLACEMENT 
GAGE  ASSEMBLY  AND  AMPLIFIER  GAGE  POSITION  POINTER 


Fig.  14c  DIFFERENTIAL 
TRANSFORMER  PROBE 


Fig.  15  SURFACE  DISPLACEMENT  GAGES 
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Step  2 


With  the  tunnel  model  in  position,  soil  was  placed  to  a  height  of 
one  tunnel  diameter  above  the  crown  (  12-in.  lift)  and  compacted  by  penetration 
of  the  probe,  to  a  depth  of  one  tunnel  diameter  below  the  invert  18-in» 
penetration  in  the  above  pattern  exclusive  of  positions  1,  2  and  3. 

Step  3 

Soil  was  placed  to  containment  vessel  flange  (9-in.  lift)  and 
compacted  by  penetration  of  probe  to  tunnel  crown  in  the  pattern  shown  in 
step  1. 

Density  measurements  made  after  each  step  indicated  that  the 
placement  technique  gave  uniform  densities. 

C.  Test  Results 

A  total  of  six  tests  were  performed  on  the  two  steel  and  single 
aluminum  tunnel  model  (two  tests  per  model).  Table  3  lists  the  tunnel  model 
properties  and  the  average  soil  density  for  each  test.  After  preparation  and 
placement  of  each  tunnel  model  in  the  containment  vessel,  as  described  in  sub¬ 
section  III  B,  the  first  test  (denoted  bv  the  A  suffix)  was  executed  by 
pressurizing  the  sand  surface  in  10  and  20  psi  increments  up  to  the  maximum 
pressure  of  100  psi.  Readings  from  all  instrumentation  were  recorded  for 
each  pressure  level  during  both  the  loading  and  unloading  cycle.  Approximately 
one  day  after  the  first  test,  without  disturbing  or  recompacting  the  sand,  the 
procedure  was  repeated  to  yield  the  second  test  results  (  B  suffix)  .  Thus, 
two  tests  were  performed  on  each  tunnel  model.  Applied  surface  pressure- 
displacement  curves  are  presented  for  each  dial  gage  for  test  IB  in  figure  17. 
Gages  1  through  4  are  surface  displacement  gages  (shown  in  figure  15),  located 
zero,  2,  6  and  10  in.  from  the  center  and  normal  to  the  tunnel  center  line. 
Gages  5  and  6  give  the  vertical  tunnel  motions  at  the  points  of  protrusion. 

The  surface  displacement  variation  indicated  by  the  curves  (gages 
1  through  4,  figure  17)  is  qualitatively  representative  for  all  tests  conducted 
(  i.  e. ,  the  largest  displacements  occurc-d  at  the  center  and  outer  gage 
points)  ;  however,  the  largest  percentage  difference, between  maximum  and 
minimum  displacement, occured  for  Test  IB. 
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*  Dry  Ottawa  Sand 
+  6061  -  6T,  Yield  Stress  =  35  ksi 


Tunnel  Displacements  Surface  Displacements 


STEEL  TUNNEL  MODEL, 


Curves  of  pressure  versus  average  surface  displacement  (average 
of  Gages  1  through  4)  are  given  in  figure  18  for  each  test.  The  average  vertical 
tunnel  model  displacement  (  average  of  Gages  5  and  6)  is  also  shown  in  the 
figure  18  for  each  test.  It  is  possible  to  estimate  the  tunnel  axis  curvature 
from  radial  displacement  gage  measurements  when  the  probe  is  vertically 
oriented.  These  measurements  yield  the  middle  ordinate  deflection  of  the 
axis  fro  a  straight  line  connecting  the  ends  of  the  tunnel.  This  ordinate 
plus  the  average  displacement  of  Gages  5  and  6  represent  the  rigid  body 
vertical  displacement  component  of  the  center  section  (Section  2,  figure  11)  . 

A  maximum  middle  ordinate  deflection  of  0.012  in.  was  observed  at  the 
100-psi  pressure  level  for  Test  3A.  This  displacement  is  20  percent  of  the 
average  end  displacement  indicating  that  the  rigid  body  displacement  at  the 
center  section  is  slightly  larger  than  the  displacement  given  by  the  average 
of  Gages  5  and  6. 

Diameter  changes  from  the  initial  position  (zero  surface  pressure) 
for  each  test  are  given  in  tables  4a  through  4f.  These  values  were  found  from 
radial  displacement  gage  measurements  taken  at  15  deg  intervals  for  each 
surface  pressure  level.  The  total  diameter  change  for  each  position  was 
found  by  adding  the  set  of  measurements  obtained  180  deg  apart.  Figure  19 
shows  average  radial  displacements  for  three  surface  pressure  levels 
(pressure  increasing)  for  Test  IB. 

Bending  and  direct  stress  at  the  crown,  45  deg  below  crown, 
springing  line  and  the  invert  are  presented  for  each  test  in  tables  5a  through 
5f.  These  stresses  were  computed  from  the  strain  measurements  taken  from 
the  gages  shown  in  figure  11.  The  hoop  stresses  presented  in  the  table  have 
been  corrected  to  eliminate  the  longitudinal  strain  effects  from  the  variation 
of  longitudinal  loading.  These  variations  arise  primarily  from  the  containment 
vessel  wall  friction.  Thus  the  listed  stresses  are  those  resulting  from  loads 
that  do  not  vary  along  the  tunnel  axis. 

Consolidation  tests  were  run  on  the  Ottawa  sand  used  in  the 
experiments.  The  objective  of  the  tests  was  to  obtain  appropriate  values  of 
the  confined  modulus  (  E£)  of  the  soil.  Two  load  cycles  with  peak  stress  of 
100  psi  were  applied  to  the  soil  sample.  The  soil  density  was  109  pcf.  The 
following  ratios  of  stress  to  strain  (Ec)  during  loading  were  obtained. 
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Average  Tunnel  Displacement  Ave rage  Surface  Displacement 

(Gages  5  and  6)  (uages  1,  2,  3  and  4) 
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19  RADIAL  DISPLACEMENTS  FOR  1  5--DEG  SEGMENTS 


Table 


41 


o 

^H 

X 

co 

W 

X 

o 

z 

I-H 

H 

£ 

w 

2 

w 

o 

c 

J 

cu 

CO 

HH 

Q 

U 

i— < 

trj 

H 

W 

2 

< 


o 

vO 

ii 

X 

n) 


(U 

TJ 

O 

2 

% 

C 

C 

3 

H 


0) 

<u 


a  w 


<D 

H 


no 

00 

m 

O' 

m 

m 

m 

i-H 

O 

o 

in 

m 

no 

l 

m 

i 

o 

i 

xO 

1 

i 

H 

1 

no 

i-H 

i-H 

oo 

m 

ro 

m 

m 

i-H 

O' 

sO 

O 

m 

i-H 

O' 

vO 

m 

O' 

00 

m 

M 

c 

^H 

i-H 

i-H 

1 

I 

i-H 

1 

H 

1 

H 

1 

i 

1 

H 

•  rH 

05 

aJ 

00 

O' 

m 

xO 

vO 

o 

O' 

m 

m 

no 

m 

<U 

o 

^H 

00 

o 

^H 

1^ 

00 

o 

00 

no 

no 

o 

*H 

(J 

m 

no 

^H 

i^ 

^H 

1^ 

no 

»— H 

1^ 

i^ 

no 

<D 

1 

1 

1 

i 

1 

1 

i 

Q 

O' 

00 

O' 

m 

00 

00 

O' 

xO 

vO 

m 

no 

o 

i-H 

i-H 

t—i 

no 

O' 

no 

o 

1^ 

m 

m 

m 

no 

no 

i^ 

1 

i^ 

i 

1^ 

1 

no 

i 

no 

i 

1^ 

1 

■ 

no 

O' 

00 

m 

00 

no 

no 

O' 

m 

m 

in 

xO 

o 

o 

no 

no 

i-H 

no 

i“H 

m 

o 

m 

fM 

no 

i-H 

^H 

no 

no 

no 

i^ 

H 

rg 

00  <D 

u  > 

°  S  g 

*  o  c 

O  Jh 

fn  Q 


Ti<  Tf  m  in  O'  O  O^OOOO^nD 
OOTfrn^HOO(\J  ^  N  N  in  N 

N  N  H  i-h  nO  N  N  H  i-h 

ill  III 


O'  O  O  no  r*-  O  rororo 

fNj  O  -H  ^  O  00  o  oo  ^ 

f\J  f\J  l-H  I-H  i-H  psj  I-H  I-H 

III  III 


m 


in  O' 

no  o 

i-h  no 


O' 


^  o 
r*-  O' 


O'  m  ^ 

O'  m 


in 


o 

O' 


o 


m 

oo 


m 


m 


o  O  vO 

^  r-H 


rO 

^H 

I 


l-H 

no  oo 


00 


vO 

m 


o  ^  vO  o  m  m  rosooc 

^  oo  ^  m  r-  o'  O'  m 

ii  iii 


vO 


m 

m 


O' 

O' 


m  m  m  co  no  no 
m  m  m  m  in 

_ i  i  i 


no  no  oo  <n  O'  o 

in  in  m  no 

i  i  i  i _ 


u 

<1) 

>  2 


o  m  o 

m  r*<  vO 


m 


u 

o 

X 

o 

O' 


m 

o 


o 

no 


m  o 
m  m 


m 

nO 


u 

<D 

•M 

<u 

6 

flJ 

■  rH 

T) 

c 

•  rH 

<D 

02 

<* 

a; 

u 

u 

<D 

<D 

•*-> 

<t3 

O 

•rH 

T) 

.3 


o 

Z 


I 

I 

I 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


-  70  - 


I 

I 

I 

I 


1 

I 


I 

I 

I 

I 


if 


X 

w 

W 

X 

o 

X 

HH 


$ 


43 

«J 

H 


H 

X 

W 

S 

w 

o 

< 

pj 

w 

HH 

Q 

U 

i— i 

o$ 

H 

w 

2 

S 

Q 


armour 


m 

m 

cO 

O' 

O' 

rg 

00 

m 

v£> 

rg 

o 

co 

cO 

H 

, 

1 

rg 

cO 

l 

cO 

l 

cO 

i 

cO 

co 

00 

if 

rg 

rg 

O' 

00 

Is- 

cO 

m 

00 

o 

m 

CO 

m 

i-H 

rg 

cO 

eg 

vO 

i-H 

H 

H 

i-H 

^H 

— i 

pq 

oo 

1 

it 

1 

1 

1 

1 

I 

c 

•  rH 

CO 

o 

if 

go 

m 

o 

cO 

cO 

<TJ 

o 

o 

oo 

i-H 

o 

m 

vD 

go 

O' 

o 

CD 

cO 

pg 

H 

^H 

i-H 

i-H 

^H 

^H 

1 

i 

1 

1 

1 

1 

I 

- 

<D 

p 

sD 

cO 

vD 

r*- 

if 

o 

o 

O' 

cO 

O 

rg 

o 

O' 

i-H 

o 

O' 

o 

o 

i-H 

o 

m 

rg 

rg 

, 

^H 

1 

^H 

1 

^H 

1 

rg 

, 

i-H 

rg 

O' 

o 

o 

00 

o 

go 

O' 

o 

o 

o 

•  iH 

o 

o 

o 

rg 

^H 

oo 

o 

oo 

o 

CO 

rg 

CO 

rg 

rg 

■— i 

i-H 

i-H 

rg 

^H 

^H 

rg 

o 

a 

i 

1 

1 

1 

1 

vO 

II 

aT 

— 

— — 

43 

3 

(fl 

CD 

00 

o 

o 

rg 

oo 

eg 

if 

cO 

m 

rg 

CO 

m 

rg 

^H 

i-H 

rg 

O' 

i-H 

r“4 

o 

co 

0) 

I— i 

rg 

rg 

pH 

i-H 

•— i 

rg 

rg 

l“< 

H 

rg 

<D 

Jh 

fc 

1 

1 

i 

1 

i 

1 

l 

TJ 

O 

o 

o 

2 

co 

O' 

rg 

if 

vO 

'O 

cO 

O' 

<3 

O 

rg 

O' 

O' 

»-H 

o 

n- 

oo 

o 

^H 

O' 

I-H 

00 

rg 

i-H 

i-H 

i-H 

^H 

i-H 

i-H 

^H 

1 

| 

1 

1 

1 

(3 

a 

c 

d 

cc 

H 

if 

^H 

if 

o 

•<* 

^H 

o 

^H 

cO 

o 

O' 

r- 

00 

^H 

O' 

^H 

oo 

oo 

o 

r- 

»-* 

>sO 

^H 

<u 

0£ 

<L 

c 

•H 

cc 

CO 

flj 

If 

vD 

O' 

o 

oo 

oo 

O 

oo 

rg 

ffl 

^H 

o 

40 

m 

cO 

H 

sD 

^H 

rg 

^H 

^H 

^H 

^H 

4-* 

c 

CD 

V 

H 

vO 

if 

O' 

rg 

o 

CO 

in 

O 

go 

O' 

rg 

o 

O' 

if 

m 

00 

O' 

O' 

in 

o 

^H 

1 

i 

t 

1 

1 

i 

rg 

o 

o 

rg 

00 

oo 

^H 

O' 

vO 

rg 

If 

^H 

''O 

rg 

1 

cO 

l 

in 

in 

i 

go 

i 

rg 

i 

cO 

SO 

^  'u 

w> 

<u 

Q 

1)  Q) 
>  *» 

<u 

S  £ 

4-» 

u 

N 

•H 

u 

0 

h  <rt 

0) 

in 

o 

m 

o 

m 

X 

in 

o 

m 

o 

in 

^£2 

> 

cO 

if 

vO 

o 

o 

^H 

rg 

^H 

cO 

^H 

in 

^H 

xO 

^H 

w' 

o 

RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


-  71  - 


Note;  +  Indicates  decrease  in  diameter. 


Table  4c 
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Note:  +  Indicates  decrease  in  diameter. 
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Note:  +  Indicates  decrease  in  diameter. 


Table  5a 


DIRECT  AND  BENDING  STRESSES,  PSI 
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DIRECT  AND  BENDING  STRESSES,  FSI 


+  indicates  compression  on  outside  fibers. 


Table  5e 


indicates  tension 


+  indicates  tension 


Stress,  E  (  1st  Cycle)  E  (  2nd  Cycle) 

(  psi)  C  (  psi) _  C(  Pgi) _ 


10 

3700 

4450 

20 

3750 

5050 

40 

4500 

5850 

60 

5450 

6700 

80 

6400 

7600 

100 

7500 

8800 

D.  Discussion  of  Results 

According  to  the  theory  presented  in  Section  II  H,  a  knowledge  of 
the  apparent  modulus  (  E*  )  ,  foundation  modulus  (  k)  ,  biaxiality  ratio  (  )  , 

and  the  stress  (O^)  together  with  the  shell  properties  (E  and  a/h)  will  enable 
one  to  determine  the  stresses  in  and  the  displacements  of  a  buried  cylinder. 
The  value  of  CT^  at  the  point  where  the  tunnel  radial  displacements  and  stresses 
were  measured  is  assumed  equal  to  the  surface  pressure  for  all  calculations 
in  this  section.  Direct  determination  of  0^  and  [2>  requires  the  measurement 
of  lateral  and  vertical  soil  pressures;  however,  at  present  no  acceptable  gage 
exists  for  measuring  static  soil  pressure. 


ka 

From  equation  57  it  is  noted  that,  — — 

Eh 


^  9000,160)  =  o.054 
10 

ka? 


when  ka  ^  E  ,  for  all  of  the  tunnel  models..  Thus  neglecting  — —  in  equation 
57  and  solving  for  the  direct  stress  at  the  tunnel  crown  (  0  =0°f'1and  springing 
(0  =  90°)  and  then  adding  and  subtracting  the  resulting  two  equations  give, 
f D  (  0  ° )  +  f D  (90°)  ka 

U7K -  =  1  +  P  +  I*  ’  (62a) 


°"o 


and 


fD  (°°)  -  2  fD  (90°) 


,  ka 

“  lTF* 


ka 

E* 


05 


(  62b) 


where  f^  (0°)  is  the  crown  direct  stress.  From  these  equations,  it  is  possible 

/^L, 

E  and  since  all  other  values  are  known  for  the  test  series 

(assuming  C T  equals  the  surface  pressure).  Presumably,  the  values  of  ka/E* 
and  /3  are  constant  for  a  given  surface  pressure  level  for  each  tunnel  model, 
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however,  attempts  to  evaluate  the  quantities  from  equations  62  result  in 
rather  wide  variations.  The  discrepancies  occur  primarily  because  the 
soil  is  not  in  a  state  of  uniaxial  strain.  This  is  evident  since  the  rigid  body 
tunnel  motions  (Figure  18)  are  only  about  1  / 3  the  surface  displacements 
instead  of  1  / Z  (  Tunnel  located  at  mid-depth)  .  Also,  the  average  stress  to 
strain  ratio  of  the  soil,  based  on  surface  displacements  and  pressures, 

A 

results  in  values  of  E  that  are  2  to  4  times  those  measured  in  the  con- 
solidometer.  The  conclusion  is  that  there  exists  appreciable  shear 
stresses  at  the  soil-containment  vessel  interface. 

The  primary  source  of  error  encountered  in  attempting  to  solve 
for  and  ka/E*  from  equations  62,  arises  from  the  evaluation  of  the 
parameter  E^x^/o^.  This  parameter  represents  (2nd  mode  tunnel  displace¬ 
ment  ratio/iree- field  strain) which  cannot  be  obtained  with  sufficient  precision 
since  the  free-field  strain  at  the  tunnel  is  not  accurately  known. 

From  equation  62a  we  observe  that 

fD(0°)  +fD(9°°> 

R  crQ'a/h 

should  be  nearly  the  same  for  each  test  if  and  ka/E  are  invariant.  From 
the  measured  values  of  direct  stress(table  5)  at  the  20,  60  and  100-psi 
surface  pressure  levels,  the  following  values  for  the  ratio  (  R)  are  obtained: 


Test 

Surface  F 

ressure 

20  psi 

60  psi 

iA 

1.  56 

mm 

mm 

1.  52 

IB 

1.  56 

mm 

1.  48 

2A 

i.  54 

1.  34 

1.  31 

1.  40 

2B 

1.41 

1.  22 

1.  21 

1.  28 

3A 

1.  26 

1.  23 

1.  22 

1.  24 

3B 

1.  34 

i.  23 

1.  20 

1.  26 

Overall  Average  (R)  =  1.36 
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All  values  are  within  15  percent  of  the  overall  average  value.  The  accuracy 

of  the  value  R  is  somewhat  questionable  in  view  of  the  assumption  that  the 

vertical  free-field  pressure  (  0^)  equals  the  surface  pressure.  It  is  also 

noted  that  both  the  ratio  R  and  the  average  soil  surface  deflections 

(see  Fig.  18)  decrease  with  progressive  tests.  This  indicates  that  the 

free-field  vertical  pressure  ft"  ,  at  the  tunnel  model  instrumentation  section, 

*  o 

decreases  with  progressive  tests  due  to  increased  friction  at  the  contain¬ 
ment  vessel  wall.  Wall  friction  is  obviously  a  function  of  the  Teflon  lining 
condition  which  had  increased  wrinkling  with  each  test.  Based  on  these 
observations  it  seems  plausible  that  the  value  of  0~Q,  used  in  calculating 
the  ratio  R,  becomes  increasingly  smaller  with  each  test.  This  effect  would 
yield  a  more  nearly  constant  value  for  the  ratio  R  than  calculated  above. 
Thus,  it  appears  that  this  aspect  of  the  theory  (  R  =  constant  for  a  given  soil 
is  acceptable. 


According  to  equation  57,  the  direct  stress  at  the  crown  and  invert 
(0=0,  and  180  deg)  are  equal.  Examination  of  table  5  reveals  that  the 
crown  stress  is  larger  for  Tests  1  and  2  and  the  invert  stress  larger  for 
Test  3.  This  is  probably  due,  in  part,  to  the  manner  in  which  the  stresses 
were  corrected  to  eliminate  the  longitudinal  loading  effects.  Generally  the 
differences  are  less  than  15  percent  at  the  100-psi  pressure  level.  It  is  also 
noted  that  the  vertical  diameter  shortening  is  larger  than  the  horizontal 
increase.  Some  difference  is  expected  from  theory  since  the  radial  displace¬ 
ment  is  given  by 


w  = 


a  4-  a,  cos  2  0 
o  2 


where 


a  — 
o 


°o  a 


2  Eh 


(n-  /3+  -*2.) 


Thus  a  difference  in  orthogonal  diameter  change  of 


4a 

o 


theory.  For  =  100  psi  and  assuming  1  4-  /3 
mum  diametric  displacement  differences  (  4aQ 


can  be  accounted  for  from  the 
ICcL 

4-  Z  the  following  maxi- 

)  for  each  test  are  reasonable. 
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Test 


max  4a  (in.  ) 


1A,  B  (Steel  model,  a/h  =  60) 

2A,  B  (Steel  model,  a/h=  40) 

3A,  B  (Aluminum  model,  a /h  =  60) 


24  x  10" 
16  x  10" 
72  x  10" 


The  value  72  x  10  is  not  exceeded  for  the  most  flexible  tunnel  model 

-4 

(tests  3A  and  B)  when  the  surface  pressure  equals  100  psi.  The  value  24x10  is 
slightly  exceeded  for  the  "near"  vertical-horizontal  diameters  for  tests 
1A  and  B.  The  largest  variations  from  the  value  were  obtained  in  the  tests 
on  the  model  with  the  highest  bending  rigidity  (  Test  2A  and  2B)  . 

It  is  interesting  to  note  that  the  Iowa  Formula  (11)  for  flexible  pipe 
culverts  with  a  bedding  angle  of  180  deg  gives 


x2  - 


l.  5  <r 


18  E  I 


+  1 .  1  ka 


as  compared  to  equation  60 

(T0(l-  /3  +ka/E*) 


18  El 


+  2  ka 


The  expression  for  ka/E*  from  equation  60  is 


-  /9 


1  8  E  I 


(62c) 


Assuming  a  free-field  strain  (  cr  / E*)  at  the  tunnel  instrumentation  section 
equal  to  the  average  strain,  based  on  the  soil  surface  displacements,  there 
results, 

i .  a 

ka  '  T3  0-o 

t4  =  - r-r-r - r— 1 — 2 -  •  <62d> 

E*  M  ~  ^90 

A 
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where 


^  =  vertical  tunnel  diametric  deflection  (  Table  4,  0  =  0  deg) 

&90  =  horizontal  tunnel  diametric  deflection  (  Table  4,  9  -  90  deg) 
Asur  =  average  surface  displacement  (Fig.  18) 

The  third  term  in  the  nhmerator  of  equations  62c  and  62d  is  small  compared 
to  unity  for  all  tests  and  can  be  neglected.  The  denominator  has  the  following 
maximum  and  minimum  values  for  increasing  and  decreasing  surface 
pressures  greater  than  20  psi. 


Test 

6(  6q  - 

O 

O' 

VX> 

^sur  | 

Max 

Min 

1A,  B 

.  73 

.46 

2A,  B 

.91 

.  73 

3  A,  B 

1.45 

.  74 

The  average  denominator  value  for  all  tests,  with  pressure  greater  than 
10  psi,  is  0.854.  Since  ka  ,  E*  and  ^/3  are  soil  properties  and  the  numerator 
of  equation  62d  is  ^  1  -  ,  the  value  6  (  cbQ  -  ^90^/Asur  anc^  ^ence 

the  numerator  if  equation  62d  should  remain  invariant  for  the  test  series, 
according  to  the  theory.  It  is  believed  that  the  variation  in  the  denominator 
(  indicated  by  the  tabulated  maximum  and  minimum  values)  is  the  direct 
result  of  the  assumption  made  for  the  free-field  strain  (  G^/e*)  .  In  reality, 
the  free-field  strain  at  the  tunnel  instrumentation  section  is  certainly  larger 
than  the  average  strain  and  the  ratio  of  free-field  strain  at  the  tunnel  to 
average  strain  probably  increases  with  increasing  wall  friction.  If  these 
hypotheses  are  assumed,  the  numerator  of  equation  62d  would  tend  to  a 
more  stable  value. 
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SYMBOL  INDEX 


Section  III 


a 

a 

n 

E 

E 

E* 

E 

c 

fD<°> 

fB«» 

h 

k 

w 

X 

n 

ft 

9 

cr 


shell  radius 

peak  shell  modal  displacement 
Young's  modulus  of  shell 
E  /( 1  -\)Z) 

apparent  modulus  of  soil 

soil  modulus  in  confined  compression 

shell  static  direct  stress 

shell  static  bending  stress 

shell  thickness 

foundation  modulus 

shell  radial  displacement 

a  /a 
n 

free -field  stress  biaxial  ratio 

angular  coordinate  from  crown 

peak  free -field  pressure  in  x-direction 


IV. 


EFFECTS  OF  DIFFERENTIAL  HORIZONTAL  SOIL 


DISPLACEMENTS 

A.  Introduction 

The  analysis  given  in  this  section  can  be  used  to  estimate  long¬ 
itudinal  bending  and  shear  stresses  in  silo  linings  resulting  from  given 
differential  horizontal  soil  displacements.  The  rigidity  of  the  silo  lining 
lies  between  two  obvious  extremes,  infinite  and  zero  rigidity.  Alining  with 
no  rigidity  will  conform  to  the  free-field  displacements  as  shown  in  figure  20b 
and  no  longitudinal  bending  or  shear  stresses  would  be  introduced  while  a 
lining  with  infinite  bending  rigidity  remains  straight  and  does  not  conform 
to  the  free-field  displacements.  The  latter  approach  clearly  yields 
maximum  stresses.  The  following  investigation  uses  the  energy  method  and 
treats  the  general  case  of  silo  linings  that  possess  both  bending  and  shear 
flexibility.  This  analysis  follows  closely  the  method  presented  in  reference  2. 

The  primary  addition  is  the  inclusion  of  rigid  body  shell  rotation  to  account  for  an 
unsymmetrical  free-field  displacement  profile. 

B.  Analysis 

If  a  lining  has  longitudinal  rigidity  it  will  tend  to  remain  straight 
and  not  conform  to  the  free-field  profile.  It  is  assumed  that  the  force  exerted 
by  the  soil  on  the  lining  is  directly  proportional  to  the  relative  radial  displace¬ 
ment  of  the  lining  and  the  corresponding  free-field  particles  (  elastic  foundation)  . 

The  displacement  of  the  lining  is  obtained  by  assuming  expressions 
for  orthogonal  displacement  components  that  contain  undetermined  coefficients 
that  satisfy  the  boundary  conditions.  The  free-field  deviations  from  a  straight 
line,  AB,  (Fig.  20c)  can  be  represented  as  a  Fourier  series.  The  shell 
displacement  coefficients  are  then  obtained  by  minimizing  the  strain  energy  of 
the  shell-foundation  system. 

The  radial,  tangential  and  axial  displacements  of  the  middle  surface 
of  the  shell  shown  in  figure  21  are  assumed  as 
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"VvYli 


F roc  -  I  ic Id 
Profi  li- 


b.  Deformed  S  ilo  with 
no  Bunding  Rigidity 


;  p]  a  refl  to  f  Si  I  o 


■i  4  ri  v  1  Y"  4  Hill 

f*  1  L.  ri  sin  — 

°  m  L 


Silo  Cc* nt c-  r  Line*  Profiles 


E VATIONS  OF  SILO  LINING 


o© 


w  =  (ao  +  OCx)  cos  0 


.  1 


m,  n  -  1 
o© 


v  =  (a  +  oc  x)  sin  0  + 
o 


u  =  oc  a  cos  0  + 


o© 

I 


m  Tf  x 

a  sin  — T -  cos  n  0 

mn  L 


m,  n  =  1 


mn  m  77  x 

-  sin  — T -  sin 

n  L 


n  0  > 


m,  n  =  1 


v  m  7T  x  „ 

b  cos  - . cos  n  0 

mn  L 


(63) 


where  oC  ,  a  ,  b  and  a  are  undetermined  coefficients, 
o  mn  mn 

The  undisturbed  radial  displacements  of  the  free -field  at  a  position 
corresponding  to  the  middle  surface  of  the  shell  is 

o© 


I 


.  a  m  JT  x  a 

w  =  /  A  sin  - t —  cos  0 

o  L — ■ .  m  L 

m  =  1 


(64) 


where  are  given  free -field  displacement  amplitudes  from  the  straight 

line  AB  connecting  the  shell  end  sections  (figures  20b  and  20c). 

An  expression  for  the  strain  energy  of  the  shell  is  given  (10)  as 

I . V  ,  1  -V 


it  _  E  ha 
s  2 


r  L  rZTT 

1 

2  , 

(ve-'")2  .  ZV 

u  + 

- 7 -  +  — 

"6 

-A. 

U  CL 

a 

a  rl 

j  rzT 

r 

E  A  f 

f 

■  2  (wee+^) 

J 

o 

J 

o 

W  +  - 1 - 

xx  4 

l  a 

u  2 


dxd0  +- 


zV  .  , 

+  T  w*x  (wee+w>  + 

a 


(65a) 


*  ^  -f 

a 


^>2 

2a  ' 


dxd0 


where 


u  = 


du 


Sx  ’  Wx0  -  dx  6 


,2 
d  w 


,  etc. 
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Substituting  the  expressions  for  the  shell  displacements  (equations 
63)  into  equation  65a  and  performing  the  indicated  operations  gives  the 
following  shell  strain  energy  expression, 


U  = 
s 


EhaLTT 


'  Z_  ZLld  b2  +  J_-ld(77m  a  _nb  2  L_ 
L —  ,  T  2  mn  2  nL  mn  a  mn  J 

m,  n=  1  u  L 


g.h4---  Z 

m,  n=  1  LL  L  a  a  L 


(65b) 


x  C™  7T  mnTCx  n  ,  1  2 

2  2nL  L  amn  2  a  mn 
a 


The  strain  energy  of  the  foundation  is 

rL  r2.TC 


_  ka 
UF  ~  ~T 


r  n 

o  o 


(w-w  )  dxd9, 
o 


(66a) 


or,  upon  substitution  of  equation  63  and  64, 


8  a  +  4  cx.  L, 


TT  ka  L7T  -,/oc.L  ,  ,  T  ,  2Vl  ao  + 
UF=  — 5 2  ( 3—  +  <X  a0  L  +  a  0)  +  ff 


a  ,  -  A 
ml  m 


m=  1,  3,  5 


(66b) 


4aL 


a  ,  -A 
ml  m 


m=  2,  4,  6 


m,  n=  1 


(a2  -2a  , A  -A2  ) 

mn  ml  m  m 


The  total  energy  is 


U  =  U  +  U-. 
s  F 


(67a) 
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The  undetermined  shell  displacement  amplitudes  are  evaluated  by 
application  of  Castigliano's  Theorem  of  Least  Work,  thus 


dU 


dU 


du 


doc  b a. 


da. 


mn 


dU 


=  0  m,  n  =  1,2 


(67b) 


mn 


From  equations  65b,  66b,  67a  and  the  last  two  equations  of  67b  we 

obtain  a  set  of  linear  algebraic  equations  in  a  and  b  .  For  n  >  1 

°  mn  mn 

the  equations  are  homogeneous  and  the  result  is 


a  =b  =0,  n  1  , 

mn  mn 


(68) 


For  n  =  1  and  dropping  the  n  subscript  since  it  is  always  unity, 
we  obtain,  after  algebraic  manipulation  to  eliminate  b  , 


m^77a  h2  .  m^7 ah2  .  2  kaL 

l  1  +  - -)  +  - , -  + 


12a2'  T  6(1  -V)T^ 


1  + 

m  7Ta  (1  -V)\  12a 


h2  N 


TT e  h(i  -v) 

\ 

2) 


2L 


00 


'a  + 

m 


2  kaL 


m 


TT  n  =  1,  3,  5 


12a 


a  -  A 
n  n 

n 


>  a 


m 


m  Odd 


TTe  h  (1-y) 


(69a) 


A  + 
v  m 


4L 


m 


TT 


Z  fs 


E  ha  ( 1-1/)  ,,  h“  .  L 
kaL -  (1+  - r)  +7T 


12a 


OO 

V 

2-  _ 

n=  2,  4,  6  n 


a  -A 
n  n 


,  m  Even 


(69b) 
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Equations  69a  and  69b  respectively  produce  two  sets  of  simultaneous 
equations  sufficient  to  determine  all  a^  (m  =  1,  2,  3  .  .  .  )•  It  is  important 
to  note  that  if  all  even  (or  odd)  Am  are  aero,  all  even  (or  odd)  am  equal 
zero. 


From  the  first  and  second  of  equations  67b  and  with  all  am 
determined  from  equations  69a  and  69b  we  find  the  rigid  body  coefficients 
of  the  shell  displacements, 


co 


(69c) 


(69d) 


Also 


b  = 
m 


d+  — ^ 
I2a^ 
71  1  2"" 


-)  m7Ta  (1  -V) 


+  (1  -y)  L  (I  +  h‘ 


m 


m  =  1,  2, 


12a 


(69e) 


The  shell  displacements  can  now  be  obtained  from  equations  63  where  sum¬ 
mation  over  n  is  neglected  or  n  is  set  equal  to  unity.  The  resulting 
expressions  for  displacement  can  be  used  to  determine  stresses  directly; 
however,  several  displacement  modes  are  generally  required  to  accurately 
determine  stresses.  Obviously  the  solution  of  the  simultaneous  equations 
(69a  and  69b)  will  require  a  major  effort  if  a  large  number  of  coefficients 
(a^)  are  required.  The  work  required  to  find  stresses  accurately  is 
significantly  reduced  by  employing  the  computed  shell  displacements  to 
determine  the  shell  loading  using  only  a  limited  number  of  the  lower  dis¬ 
placement  modes.  The  shear  and  moment  are  then  readily  determined 
from  this  loading. 
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The  radial  pressure  acting  on  the  surface  of  the  shell  is 


P 

P 


k  (wq  -  w) 


0-0 

I 


n  =  1 


(A  -a  )  sin 
'  n  n' 


n7Tx 

~U~ 


cos  6. 


The  resultant  lateral  load  per  unit  length  of  the  shell  is 


F 


cos  0  d  0 


F  =  77ke 


oo 

I 


OCX  -  a  +  / _  (A  -  a  )  sin  — T - 

o  1 — •  ,  n  n'  E 

n  =  1 


The  shear  and  bending  moment  are  found  from 


-f 


V  =  |  F  dx 
o 


V  =  7Tka  —oc  — -  a  x  + 


oo 

L  /  An  an  , ,  n  7Tx\ 
^  ’  77  nT  1  ”  1  "  °0S 


and 


'•x 


M  =  I  V  dx 


M  =7Tka 


aoT 


OO 

2  .  )  A  -a 

*  +J=L  . -2L  -n(x  - 

2  7T  n  =  1  n 


L  n  77x  x 
sm  — ▼ — ) 


n7T 
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The  elastic  extreme  fiber  bending  stresses  and  the  maximum  shear 
stress  for  a  reinforced  concrete  silo  lining  with  no  axial  load  can  be  found 
from  the  following  formulas:  (Uniform  distribution  of  the  reinforcing 
steel  at  each  face  and  a  distance  of  4-in.  between  the  center  of  the  bar  and 
the  concrete  surface  is  assumed) 

Me 

Q 

f  =  — T —  ,  Concrete  Stress 

c  I 

nMc 

g 

f  =  — j - ,  Reinforcing  Steel  Stress  (71) 


T  =  k  ,  Shear  Stress  (Cracked  Section) 

where  (see  Fig.  22), 


c  =  a  +  -r  -  Y 


c  =  a  +  —  +  y  -  4" 

S  L. 


y  =  a  cos  0 

I  =  (0  +  sin  0  cos  0)  (a^  +  — )  ah  +7Thnp  (a"5  +  ah^  +  12  ah  + 

,2 


3  .  3  L  2 

T 


12a 

|(7T-  9)  cos  0  + 

[>  ♦  <zY  -  4] 

L  a  J 

A 

sin  0j 

n  =  ratio  of  modulus  of  elasticity  of  steel  to  that  of  concrete 
p  =  ratio  of  total  longitudinal  steel  area  to  gross  concrete  area. 


and  0  is  found  as  a  root  to  the  equation 


0  -  ( 1  +  - s  )  tan  0  +  7Tn  p  =  0 


12a 


(72) 
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A  plastic  analysis  has  been  employed  to  obtain  the  following 
expression  for  the  ultimate  resisting  moment  of  the  silo  lining.  A  constant 
concrete  compressive  stress  of  0.  85  f'  and  constant  compressive  and 
tension  steel  stress  of  f  is  assumed. 

yp 


M 


ult. 


=  2 


p  f  a^h 

yp 


(7T-9p)  cos 


e  + 

p 


l  + 


2 

(4) 


4h 

~Z 


sin  & 

P 

(73) 


where 


f  =  yield  stress  of  reinforcing  steel 

f ' 

0  =77 p/(2p  +  0.  85  -£—) 

p  yp 
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The  procedure  is  illustrated  in  the  following  example  for  a  concrete 
silo  40  ft  in  diameter,  160  ft  deep  with  a  24-in  wall  thickness.  A  homo¬ 
geneous  shell  with  E  =  3  x  10^  psi  and  \)  =  0  and  a  value  of  ka  =  5000  psi 
is  assumed.  A  typical  horizontal  displacement  profile  of  the  free -field 
is  shown  in  figure  23.  A  three  term  approximation  of  the  given  profile 
is  used  in  this  example  defined  by  the  following  coefficients, 

A^  =  0.  75  in. 

A^  =  0.  185  in. 

A  ^  =  0.113  in. 

A  =  0,  n  >3 

n 

From  equations  69a  and  69b  and  considering  only  five  displacement 
modes  we  obtain  as  the  solution  of  the  simultaneous  equations 


al 

=  0.  0860  in. 

a  .  =  0.  0158  in. 
4 

a2 

=  0.  0881  in. 

a^  =  -0.  001 1  in. 

a3 

=  -0.  0024  in. 

Substituting  these  values  into  equations  69c  and  69d, 

_  4 

OC  =  -0.888  x  10  in.  /in. 

a  =  0.5  326  in. 

o 

The  displaced  center  line  of  the  silo  lining  using  these  coefficients  is 

illustrated  in  figure  23.  The  displaced  center  line  of  a  rigid  silo  is  also 

shown.  The  rigid  silo  displacement  coefficients  (  a^  and  (X)  are  found 

from  equation  69c  and  69d  where  all  a  =0  for  n>0.  The  resultant  load, 

n 

shear,  and  bending  moment  are  shown  in  figure  2  3  for  three  flexibility 
assumptions,  (1)  the  rigid  silo,  (2)  three  displacement  modes  a^.  .  .a  ^4  0 
and  (3)  five  displacement  modes  a^.  .  .  a^  4  0.  It  is  important  to  note  that 
only  minor  differences  exist  between  the  three  and  five  mode  curves. 
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Typical  Free-Field  Profile 


Free-Field  Profile  Used  in  Example, 
(.  7  5  sinlLi.  +  _  185  sin  X  +  .113  sin 


Displaced  fc-  ,  a  +  OL 
o 


x  +  \  a  sin 

L,  n 


n  =  1 


■Displaced  of  Rigid  Silo 
(.  6777  -  .  000 1 84x) 


Fig.  23  FREE  -  FIELD  AND  SILO  HORIZONTAL  DISPLACEMENTS 
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5  Displacement  Modes 


3  Displacement  Modes 


The  ultimate  resisting  moment,  ratio  of  maximum  computed  j 

moment  to  ultimate  resisting  moment,  and  maximum  elastic  stresses  for 
three  percentages  of  longitudinal  reinforcing  steel  are  given  in  table  6  | 

for  the  example  utilizing  the  five  displacement  mode  results. 
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SYMBOL  INDEX 


Section  IV 


m 


a  ,  a 
mn  m 


b  ,  b 
mn  m 

c  ,  c 
c  s 

E 

E 

F 


f  ,f 

c  s 


yp 

h 

k 

L 

M 


P 

P 

Us’UF 
u,  v,w 

V 

w 

o 

X,  z 

oc 

e 

V 
T 


Shell  Radius 

Free -Field  Displacement  Coefficients 
Shell  Displacement  Coefficients 

Distance  from  Neutral  Axis  to  Extreme  Concrete  and  Steel  Fibers 
Young's  Modulus  of  Shell 
=  E  /  ( 1  -V) 

Lateral  Load 

Concrete  and  Steel  Stress 

Yield  Stress  of  Reinforcing  Steel 

Shell  Thickness 

Foundation  Modulus 

Shell  Length 

Bending  Moment 

Ratio  of  Steel  Area  to  Gross  Concrete  Area 

Ratio  of  Young's  Modulus  of  Steel  to  that  of  concrete 

Radial  Pressure 

Shell  and  Foundation  Strain  Energy- 

Shell  Displacements 

Shear 

Free-Field  Radial  Displacement 

Shell  Coordinate 

Shell  Displacement  Coefficient 

Shell  Coordinate 

Poisson's  Ratio 

Shear  Stress 
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V. 


AXIAL  LOADING  OF  SILO  LININGS 


A.  Introduction 

As  discussed  in  Section  1,  the  silo  lining  will  be  subjected  to  axial 
loading  generated  by  blast  pressure  applied  directly  to  the  silo  cover  and  by 
longitudinal  surface  shear  or  skin  friction  between  the  structure  and  the 
surrounding  soil  medium.  When  the  silo  cover  is  supported  directly  on  the 
silo  lining  (monolithic  cover),  axial  loads  are  transmitted  directly  to  the 
lining  by  air-blast  pressure  on  the  cover.  The  dynamic  effects  of  the  rapid 
application  of  this  cover  load  are  discussed  in  Appendix  A.  The  axial 
loading  of  silo  linings  considered  in  this  section  is  treated  as  a  static 
problem,  i.  e.  ,  the  air-blast  pressure  surface  loading  is  taken  as  a  static 
pressure.  Since  silo  linings  have  much  more  rigidity  in  the  vertical 
direction  than  the  soil  they  replace,  there  is  relative  vertical  motion  between 
the  structure  and  the  surrounding  soil  medium.  This  relative  motion, 
illustrated  schematically  in  figure  25,  generates  shearing  stresses  at  the 
silo-soil  interface  that  produce  axial  silo  loads  in  their  own  right. 

For  a  silo  with  monolithic  cover,  the  skin  friction  produces  axial 
loading  in  addition  to  that  generated  from  the  cover  loading.  When  the  silo 
cover  is  supported  on  a  foundation  independent  of  the  silo  lining,  the  skin 
friction  is  the  only  source  of  axial  loading.  Often  the  axial  loads  are  suf¬ 
ficiently  large  to  control  the  thickness  of  the  lining.  An  emperical  method 
that  is  commonly  employed  for  determining  skin  friction  magnitudes 
involves  two  simplifying  assumptions;  1)  the  normal  stress  acting  on  the 
lining  is  the  same  as  the  free-field  horizontal  stress,  2)  the  friction 
stress  is  found  by  multiplying  the  normal  stress  by  the  tangent  of  the 
soil's  angle  of  internal  friction.  The  skin  friction  analysis  presented  in 
this  section  considers  the  equilibrium  of  a  conical  wedge  of  soil  adjacent 
to  the  silo  that  is  in  a  state  of  impending  downward  motion  relative  to  the 
silo.  Expressions  are  derived  for  upper  bound  values  of  the  surface  shear 
stress  for  silo  structures  with  both  independent  and  integral  covers. 

Push-out  tests  of  walls  embedded  in  sand  were  presented  in  the 
first  interim  report  (  7)  to  determine  the  usefulness  of  the  upper  bound  value 

ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


107  - 


Ground  Surface 


given  by  the  theory.  Good  comparison  was  obtained  between  test  results 
and  the  theory.  Similar  push-out  tests  have  been  performed  on  a  small  silo 
model.  Comparison  between  the  measured  and  theoretical  skin  friction  is 
fair  for  the  silo  model  push-out  test  series  but  cannot  be  considered  as  entirely 
satisfactory  in  view  of  the  difference  encountered;  i.  e.  ,  the  measured  friction 
force  was  approximately  20  percent  higher  for  some  of  the  tests  than  the 
value  given  by  the  upper  bound  theory.  Tests  on  a  segmented  silo  model 
with  surface  pressure  up  to  275  psi  are  also  given.  These  tests  were 
inconclusive  since  the  axial  rigidity  of  the  model  was  too  low  for  the  soil 
employed  (  dense  Ottawa  sand)  . 

B.  Analysis 

It  is  clear  from  the  preceding  discussion  that  axial  load  will  be 
transmitted  through  the  soil  to  the  structure  by  shear,  since  the  silo  lining 
has  much  more  axial  rigidity  than  the  soil  it  replaces.  That  is,  relative 
vertical  displacement  of  the  soil  adjacent  to  the  lining  with  respect  to  the 
structure  (  see  figure  25)  causes  shear  stress  at  the  interface.  Under  loading 
there  is  a  downward  displacement  of  the  soil  with  respect  to  the  structure 
extending  from  the  surface  to  the  neutral  plane,  (the  plane  at  which  there  is 
no  relative  motion)  .  Below  the  neutral  plane,  the  relative  displacement 
of  the  soil  with  respect  to  the  structure  is  upward.  The  location  of  the 
neutral  plane  is  determined  by  the  requirement  that  the  structure  be  in 
equilibrium  under  the  action  of  vertical  forces.  The  axial  force  in  the  lining 
is  a  maximum  at  the  neutral  plane  and  will  be  equal  to  the  sum  of  the  loads 
transmitted  directly  by  the  cover  plus  the  sum  of  the  friction  forces  above 
this  plane.  Thus,  the  resulting  problem  becomes  the  determination  of  the 
sum  of  the  friction  forces  acting  on  the  silo  above  the  neutral  plane.  The 
following  analyses  apply  to  cases  where  the  silo  cover  is  supported 
independent  of  the  lining  (  independent  cover  case)  and  where  the  silo  cover 
is  an  integral  part  of  the  silo  lining  (  monolithic  cover  case)  . 

In  both  instances  the  analyses  involve  the  assumption  that  there 
exists  a  region  of  impending  failure  bounded  by  the  cylindrical  surface  of  the 
silo  lining  (  AB)  and  a  conical  surface  (  BC)  through  the  soil  extending  up¬ 
ward  from  the  neutral  plane  at  an  angle  9  with  the  horizontal  (  see  figure  26)  . 
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The  following  assumptions  are  also  made: 

1)  The  soil  outside  the  conical  surface  (BC)  is  in  a  state  of 
impending  downward  motion  relative  to  the  soil  inside  the 
surface.  The  motion  of  the  soil  adjacent  to  the  silo,  above 
the  neutral  plane,  is  downward  relative  to  the  lining. 

2)  The  full  shearing  strength  of  the  soil  is  mobilized  at  the  silo 
surface  (  AB)  and  the  conical  surface  (  BC)  and  is  given  by 
Coulomb's  Law, 


~C  -  c  +■  n  tan  <ji  (  74) 

where 

~C  is  the  unit  shearing  strength, 

c  is  the  unit  "no-load"  shearing  strength  or  unit  cohesion, 
n  is  the  applied  normal  stress  at  the  failure  surface,  and 

$  is  the  angle  of  internal  friction  of  the  soil. 


Consider  an  elemental  wedge  cut  from  the  region  of  failure  by 
vertical  planes  as  shown  in  figure  27.  The  wedge  is  in  equilibrium  under 
the  action  of  the  following  resultant  forces. 

W  the  weight  of  the  soil  wedge  ABC 

P  the  resultant  of  the  surface  pressures  acting  on  the  face 

AC 

Fj  the  resultant  of  the  friction  stresses  acting  on  BC 

F^  the  resultant  of  the  friction  stresses  acting  on  AB 

Cj  the  resultant  of  the  cohesive  stress  on  BC 

C^  the  resultant  of  the  cohesive  stress  on  AB 

F^  the  resultant  of  the  circumferential  stresses  acting  on 
the  vertical  surfaces  of  the  wedge. 

The  surface  pressure  resultant  for  a  silo  with  an  independently 
supported  cover  is 

P  -  £  (  a  +  L  cot  6)  2  doc  •  (75) 
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This  equation  assumes  that  the  total  surface  pressure  acting  on  the  cross 
sectional  area  of  the  lining  (  elemental  area  OA,  figure  27)  is  applied  to  the 
surface  AC. 

For  a  silo  with  monolithic  cover, 

p  =  P_Lco*e -  (2a  +  L  cot  q)  doC  _  (76) 

The  soil  wedge  weight  is, 

2 

W  =  wL  cot  6  (  a  +  -  )  dot  ,  (  77) 

where  w  =  p  L,fl  and  p  is  soil  density.  The  resultant  cohesive  forces  are, 

C1  =  (a  +  4  cot  ®)  (78) 

and 


=  caL  doc  (  79) 

The  friction  forces  and  F^  can  be  specified  in  terms  of  the  un¬ 
known  force  F^,  through  the  conditions  of  vertical  and  horizontal  equilibrium. 
If  the  ratios  Fj/f^  and  are  known,  the  angle  G  that  minimizes  F^  can 

be  found.  F2  obtained  in  this  way  is  correct  provided  the  failure  surface  is 
the  proper  one.  The  determination  of  F^  requires  a  detailed  consideration 
of  equilibrium  throughout  the  wedge;  however,  an  upper  bound  on  the 
resultant  F2  can  be  found  by  assuming  F^  =  0.  For  this  case,  the  equations 
of  equilibrium  in  the  horizontal  and  vertical  directions  are 


j  cos  9  +  F  j  sin  (  0  -I-  $)  -  F2  cos  =  0 


(80) 


P  +  W  +-  CjSinG-C2-FjCos(G+f()  -  F2  sin  $  =  0.  (81) 


Eliminating  F^  from  equations  80  and  8i  and  solving  for  F2  yields, 

P  _  P+W-t-Ci  [sin  6  +  cos  6  cot  (  6  +  $)]  -  Cz 
2  ~  sin  0  ■+•  cos  $  cot  (9  +  0) 


(82) 
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The  average  shear  stress  attributed  to  friction  at  the  silo-soil  interface  is 

F?  sin  $ 

ntan*  *  IXeToT  (83> 


Substituting  equations  75,  77,  78  and  79  into  the  nondimensional  form  of  this 
equation  results  in  (  independent  cover  case)  , 

n  tan 


-J-  4-  Cot  G)  ^  4-  ^  cot  0  J  cot  (  G  4-  ft)  4- 


L 

2a 


1  4-  cot  9  cot  (  G  4-  0)|  ^  4-  —  cot  G  (  1  4-  cot^  0) 


(84) 

x 


1 


( tan  G  4-  tan  $)  sin  $  cos  ^ 

The  correct  angle  (  G)  is  that  which  produces  a  minimum  value  for  n  tan  0//p. 
The  angle  is  readily  found  as  a  root  to  the  equation, 


d  ,  n  tan  JB 

P  ' 
d  cot  0 


=  0. 


Performing  the  indicated  operation  yields 


tan  G  4-  A  tan  0  4-  B  =  0 


(85) 


where 


and 


A  =  -  — 

a 


B  =  -2 


f  + 


1  L  1  w  L  ,  /  1  ,  w  \  *.  A 

13 —  — -f-  -  ■  — —  4”  (  1  4“  — )  tan  C 
7  a  3  p  a  v  p  '  r 


[  -§■  *  ( 1  +  1 — |r  >  • 


Considerable  simplification  is  afforded  by  neglecting  the  second-order  effect 


w  w 

of  — ;  thus  when  —  =  0 
P  P 


tan  0  — 


Za 


1  +  /,  4-  8-1-  (  § 


4-  tan  fi) 


(86) 
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The  resulting  ratio  of  average  shear  stress  at  the  silo  interface  to  the 
surface  overpressure,  z /p,  for  the  independent  cover  case  is 


5  =  4-  + 


n  tan  0 
P 


Independent  Cover 


(87) 


where  the  second  term  in  the  right  hand  member  is  given  by  equation  84  with 
minor  error  involved  when  0  from  equation  86  is  used.  Curves  t/p  versus 
a/L  are  plotted  in  figure  28  for  various  values  of  $  when  w/p  =  c/p  =  0  for 
the  independent  cover  case. 

Substituting  equations  76  through  79  into  equation  83  yields,  for  the 
monolithic  cover  case, 


n  tan  0  _  1 

~ - T 


|  C 

cot  0(1+  ^  cot  0)  4-  —  cot 


0 


jcot  (  0  +  ^)  — £  1  +  cot  0  cot  (  0  +  | 


-h 


(88) 


4-  cot  0  (  1  +  cot^  0  ) 


(tan  0  +  tan  $)  sin  0  cos  jif. 


Here,  it  makes  no  sense  to  seek  the  minimum  value  of  equation  88  since  the 
smallest  realistic  value  of  n  tan  $/p  occurs  when  0=90  deg.  For  this  angle, 
the  failure  zone  can  be  thought  of  as  a  thin  wedge  with  an  infinitestimal  open¬ 
ing  angle,  in  which  case,  L1  Hospital's  rule  is  applied  to  the  indeterminate 
0:6$o  form  of  equation  88.  The  resulting  ratio  of  average  shear  stress  at 
the  silo  interface  to  the  surface  pressure  for  the  monolithic  cover  case  is 


£  + 


1  -I-  ^+-^-(^4  _  tan  0)  J  sin  2  0  ,  (  89) 


Monolithic  Cover. 


C .  Silo  Model  Push-Out  Tests 

To  determine  the  usefulness  of  the  upper-bound  skin  friction  theory 
given  in  the  previous  section,  push-out  tests  have  been  conducted  with  a  small 
silo  model  embedded  in  sand  for  comparison  of  measured  values  to  those  given 
by  equations  87  and  89.  Relative  motion  for  these  experiments  similar  to  that 
described  in  the  analysis,  i.  e.  ,  downward  motion  of  the  soil  adjacent  to  the 
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Shear  Stress  to  Surface  Pressure  Ratio,  'tr/p 


lining,  has  been  obtained  by  forcing  the  silo  model  up  through  the  sand  with  a 
hydraulic  jack.  Eight  tests  were  performed,  two  with  an  independent  cover 
and  six  with  monolithic  or  integral  cover.  Both  loose  and  dense  sand  conditions 
were  used  and  surface  pressures  from  zeroto  20  psi  were  applied  prior  to 
application  of  the  silo  model  displacement  load.  A  sketch  of  the  test  arrange¬ 
ment  is  shown  in  figure  29. 

The  tests  were  performed  in  the  3  -  ft  diameter,  3-ft  high  contain¬ 
ment  vessel  that  was  used  for  the  tunnel  model  tests  (  section  III  B)  .  The 
silo  model  was  2-ft  long,  6-in.  in  diameter  and  weighed  45  lb.  A  9-in. 
diameter  cover  was  placed  over  the  top  of  the  model  in  two  of  the  tests  to 
simulate  an  independently  supported  cover.  A  small  clearance  was  provided 
between  the  cover  and  silo  model  to  permit  relative  motion  between  the  two. 

The  remainder  of  the  tests  were  performed  without  the  independent  cover 
to  simulate  a  silo  with  a  monolithic  cover  the  same  diameter  as  the  silo. 

Twenty  V-grooves  per  inch  (  l/32-in.  deep)  were  cut  in  the  outside  face  of  the 
silo  model  normal  to  the  silo  axis  to  assure  a  value  of  $  equal  to  that  of 
sand-on-sand  at  all  failure  surfaces.  Direct  shear  tests  confirm  that  $  for 
the  grooved  surface  is  higher  than  the  sand-on-sand  value;  consequently, 
failure  occurs  on  a  cylindrical  surface,  at  the  silo  model  wall,  that  has  a 
radius  slightly  larger  than  the  model  radius. 

Dry  Ottawa  sand  with  a  grain- size  distribution  shown  in  figure  16 
was  used  for  all  tests.  Loose  sand  with  average  density  ranging  from  99  to 
100.  5  pcf  were  used  for  half  the  tests.  The  loose-sand  condition  was  obtained 
by  slowly  pouring  the  sand  in  the  pressure  vessel  with  a  hand-scoop  in  a 
uniform  manner.  Free-fall  during  placement  was  limited  to  less  than  1  inch. 
Dense  sand  was  obtained  by  penetration  of  a  one-horsepower,  10,  000  rpm 
concrete  vibrator  probe  the  depth  of  the  sand-filled  pressure  vessel  at 
several  positions  around  the  model.  Average  dense  sand  densities  varied  from 
107.  7  to  111.  0  pcf.  The  relationship  between  the  angle  of  internal  friction 
and  the  sand  density,  given  in  figure  30  was  determined  from  triaxial  test  data 
for  confining  pressures  of  10  and  x5  psi. 

Table  7  gives  the  model  arrangement,  surface  pressure,  soil 
properties  and  the  test  results  for  the  eight  tests.  Load-displacement  curves 
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Showing  Independent  Showfng  Monolithic 

Cover  Arrangement  Cover  Arrangement 


Fig.  29  SILO  MODEL  PUSH- OUT  TEST  ARRANGEMENT 


Table  7 


SILO  MODEL  PUSH- OUT  TESTS 


Test  No. 

Silo  Model 
Arrangement 

Average 

Sand  Density 
(pcf) 

0 

(deg) 

Surface 

Pressure 

(psi) 

Max 

Shear 

Load 

(lb) 

mem 

■ 

1 

Independent  Cover 

100.  5  (loose)* 

30.  0 

10 

1205 

0.  027 

0.  267 

1. 026 

a 

Independent  Cover 

1  10.  3  (dense) 

36.  5 

7.  0 

4305 

0.  220 

1. 29 

1.666 

H 

Monolithic  Cover 

99. 6  (loose) 

29.7 

10 

940 

0.  014 

0.23 

0.  491 

1 

Monolithic  Cover 

111.0  (dense) 

37.  5 

5 

3535 

0.  354 

1. 68 

1.412 

5 

Monolithic  Cover 

100  (loose) 

30.  0 

20 

1840 

0.  010 

0.  22 

0.  470 

6 

Monolithic  Cover 

108.  9  (dense) 

35.  5 

5 

3225 

0.  292 

1. 53 

1. 247 

7 

Monolithic  Cover 

99  (loose) 

29.  5 

0 

125 

- - 

Apparent  c  =.055  psi 

8 

1 

i 

Monolithic  Cover 

1 07 .  7  (dense ) 

34.  5 

0 

625 

' 

1.31  psi 

—  ose,  dense  refers  to  placement  technique  (see  text). 
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Fig.  30  ANGLE  OF  INTERNAL  FRICTION  VERSUS  DENSITY, 

OTTAWA  SAND 
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Load,  ( lb) 


Fig.  31a  LOAD  DISPLACEMENT  CURVES  FOR  SILO  MODEL 
PUSH-OUT  TESTS,  INDEPENDENT  COVER 
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Load,  ( lb) 


Fig.  31b  LOAD-DISPLACEMENT  CURVES  FOR  SILO  MODEL 
PUSH- OUT  TESTS,  MONOLITHIC  MODEL 
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for  each  test  are  given  in  figures  31a  and  3 lb.  Tests  7  and  8  (monolithic 
cover  tests)  were  run  with  no  surface  pressure  to  obtain  a  value  for  the 
"no-load"  shearing  strength  or  the  apparent  cohesion  (c)  at  the  silo  model 
surface.  The  apparent  cohesion  was  found  according  to  the  assumption  that 
the  surface  shear  ~C  given  by  the  theory  (  Eq.  89)  accurately  represents 
the  measured  average  shear  value.  Thus,  by  taking 

-nl  -  r] 

theory  measured  / 

it  is  possible  to  compute  the  value  of  c  since  p,  w,  a/L  and  0  are  known. 
The  result  is:  c  =  .  055  psi,  Test  7  (  Loose  Sand,  p  =  99  pcf);and  c  =  1.31 
psi,  Test  8  (  Dense  Sand,  p  =  107.  7  pcf)  . 

This  apparent  cohesion  is  fictitious  in  that  it  is  actually  caused  by 
intergranular  pressures  created  during  soil  placement.  These  intergranular 
pressures  generate  initial  normal  soil  stress  at  the  failure  surfaces  over  and 
above  the  normal  stress  determined  from  equilibrium  considerations.  The 
initial  normal  stresses  impart  a  "no-load"  shearing  resistance  property  that 
is  assumed  to  be  identical  to  that  created  by  cohesion.  The  apparent 
cohesive  stresses  found  from  Tests  7  and  8  are  used  to  obtain  c/p  ratios 
for  tests  1  through  6  by  assuming  a  linear  relationship  between  c  and  the 
sand  density.  The  value  of  c  /p  obtained  in  this  manner  is  given  in  table 
7  for  each  test. 

Maximum  measured  shear  loads  at  the  silo  model  interface  and 
average  theoretical  and  measured  shear  stress-to-surface  pressure  ratios 
te/p)  are  also  given  in  the  table.  The  maximum  shear  loads  were  found  by 
subtracting  the  silo  model  weight  from  the  maximum  jack  load  for  the 
independent  cover  tests  and  subtracting  the  model  weight  and  surface 
pressure  (  p)  acting  on  the  silo  model  for  the  monolithic  cover  tests. 

Measured  and  theoretical  z/p  values  for  the  loose  sand  tests 
compare  very  poorly.  The  wide  discrepancy  is  probably  caused  by  the  small 
void  ratio  reduction  resulting  from  grain  redistribution  in  the  vicinity  of  the 
silo  model  wall  which  relieves  part  of  the  normal  pressure  at  the  silo  model 

surface.  It  is  concluded  that  the  analysis  is  not  applicable  to  soils  that 
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experience  void  ratio  reductions  during  shear  (e.  g.  ,  loose  sand)  .  It  is 
interesting  to  note  that  an  average  normal  pressure  (n)on  the  model  surface 
of  l/3  the  surface  overpressure  (p)  results  in  r/p  values  of  approximately 
0.  21  which  compares  favorably  with  the  measured  ratios  for  the  loose  sand 
tests. 

For  the  dense  sand  tests,  the  x/p  values  given  by  the  theory  compare 
fairly  well  with  the  measured  values.  However,  for  the  monolithic  cover 
dense  sand  tests  (Tests  4  and  6),  the  measured  r/p  ratio  is  approximately 
20  percent  higher  than  the  value  obtained  from  the  theory.  Thus,  for  these 
tests,  the  theory  did  not  yield  an  upper-bound  value.  This  indicates  that  the 
assumed  conical  failure  surface  through  the  soil  is  incorrect. 

D.  Segmented  Silo  Model  Tests 

These  tests  simulate  a  more  realistic  silo  model  environment 
than  the  push-out  tests  described  in  the  previous  section.  For  these 
experiments,  shear  loads  developed  at  the  silo  model  surface  are  generated 
by  relative  motion  between  the  structure  and  soil.  This  relative  motion  is 
created  by  soil  compression.  The  silo  model  was  embedded  in  dense  Ottawa 
sand  contained  in  the  Foundation  Dynamic  Soil  Facility.  The  test  arrangement 
is  show  schematically  in  figures  32a  and  32b  for  the  independent  and  monolith 
cover  tests.  Two  tests  were  run,  one  with  and  one  without  the  independently 
supported  cover.  Surface  pressure  up  to  285  psi  was  employed. 

To  accentuate  the  silo  skin  friction  forces  (  generated  by  the 
relative  displacement  at  the  silo- soil  interface)  the  segmented  silo  model 
has  been  made  as  rigid  as  possible  in  the  longitudinal  direction  acknowledging 
that  the  silo  segment  loads  must  be  measurable  even  for  small  friction  forces. 
The  cantilever- Maypole  arrangement  for  the  segmented  silo  shown  in  figure 
33  was  used  for  these  tests.  The  segmented  silo  model  consists  of  four,  equal, 
6-in.  diameter  cylindrical  segments  6-in.  long,  each  independently  pin  connected 
to  three  cantilever  arms  radiating  from  a  3/4-in.  diameter  center  pole  that  is 
supported  at  the  silo  model  base  plate  (see  figure  33  and  34)  .  A  3/32-in 
clearance  was  provided  between  each  segment  and  at  the  base  and  cover.  The 
cylinder  segments  were  fabricated  from  3/16-in.  wall  steel  tubing  and-  20  V- 
grooves  per  inch,  l/3  2  -in.  deep  were  cut  in  the  outside  surfaces  to  assure  a 
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Typical 

Segment 


silo-soil  interface  angle  of  internal  friction  equal  to  the  sand-on-sand  value. 

An  assembled  and  exploded  view  of  a  typical  segment  is  shown  in  figure  35 

Silo  lining  model  shear  resultants  of  each  silo  segment  were  deter¬ 
mined  from  strain  measurements  made  with  Budd  etched-foil  gages,  type 
C6-121  and  C6-141-B.  The  location  of  the  gages  is  shown  in  figure  33.  The 
C6-121  gages  are  located  on  the  silo  segment  support  arms  near  the  center 
pole  in  order  to  register  the  highest  arm  strain.  These  gages  are  mounted 
in  pairs  on  each  arm,  one  gage  each  on  the  upper  and  lower  surfaces.  The 
C6-141-B  gages  are  mounted  in  pairs  on  opposite  sides  of  the  3/4-in.  diameter 
center  pole  below  the  silo  segment  support  arms.  These  gages  register  the 
sum  of  the  segment  loads  and  serve  as  a  check  on  the  individual  segment  loads 
obtained  from  the  arms.  The  output  of  each  pair  was  averaged  to  eliminate 
the  effects  of  local  bending  of  the  pole.  The  9-in.  diameter  independent  silo 
cover  was  the  same  as  used  for  the  push -out  tests. 

The  silo  cover  and  sand  surface  vertical  displacements  were 
measured  during  the  surface  pressurization  by  two  LVDT  s,  one  at  the 
center  of  the  cover  and  one  10-inches  from  the  pressure  vessel  wall.  Figure 
36  shows  the  two  LVDTs  in  place  with  pressure  vessel  head  removed.  The 
ARF  Dynamic  Soil  Facility,  shown  in  figure  36,  consists  of  a  pressure  vessel 
fabricated  from  a  4 -ft  ID  by  4 -ft  cylinder  with  a  fixed  dished  head  on  the 
bottom  and  a  removable  flat  head  on  top.  Design  was  based  on  300 -psi 
internal  pressure.  The  Ottawa  sand  used  for  the  tests  was  the  same  as  used 
for  the  silo  push-out  tests.  The  average  sand  density  of  108.  9  pcf,  for  both 
tests,  was  obtained  by  penetration  of  a  one -horsepower,  10,000  rpm  concrete 
vibrator  probe  in  a  symmetrical  manner  around  the  silo  model. 

Surface  pressure  versus  displacement  curves  for  the  sand  surface 
10-in.  from  the  pressure  vessel  wall  (14-in.  from  the  center  of  the  model) 
and  the  silo  cover  are  given  in  figure  37  for  both  tests.  Curves  of  pressure 
versus  shear  load  for  each  of  the  four  segments  are  given  in  figures  38a  and 
38b.  These  shear  loads  were  obtained  from  segment  arm  strain  measure¬ 
ments. 

It  is  observed  in  figures  38a  and  38b  that  the  shear  loads  are  much 
smaller  than  predicted  from  the  upper-bound  analysis  (given  in  section  V  B) 
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a.  ASSEMBLED  VIEW 


b.  EXPLODED  VIEW 


Fig.  35  TYPICAL  SILO  MODEL  SEGMENT 
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Rubber  Sealing 
Membrane 


Fig.  36a  SEALING  MEMBRANE  AND  LVDT  ARRANGEMENT 


Fig.  36b  ARF  DYNAMIC  SOIL  FACILITY 
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Fig.  37  SURFACE  PRESSURE  VERSUS  SILO  COVER  AND  SOIL  DISPLACEMENT 
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Fig.  38a  PRESSURE  VERSUS  SEGMENT  SHEAR,  SILO  WITH  MONOLITHIC  COVER 


Surface  Pressure,  (psi) 


Silo  Segment  Shear  Load,  (lbs.  ) 

Fig.  38b  SURFACE  PRESSURE  VERSUS  SEGMENT  SHEAR, INDEPENDENT  COVER 


since  the  relative  motion  of  the  soil  with  respect  to  the  silo  is  quite  small. 
This  is  apparent  from  the  monolithic  cover  curves  in  figure  37;  that  is, 
the  soil  and  cover  displacements  are  nearly  identical  whereas  the  soil 
strain  must  exceed  the  model  strain  to  obtain  the  necessary  relative  motion 
(  see  figure  24)  .  Unfortunately,  the  soil  strains  were  even  less  than 
anticipated.  The  surface  displacements,  which  are  a  measure  of  the  soil 
strains,  were  considerably  smaller  for  these  tests  than  for  the  tunnel 
model  tests  (Fig.  18).  The  untreated  wall  (no  Teflon  lining)  of  the  ARF 
dynamic  pressure  vessel  is  probably  the  primary  cause  of  the  nearly 
trivial  soil  strains.  The  conclusion  is  that  the  silo  is  more  flexible  than 
the  soil;  consequently,  the  silo  strains  match  the  soil  strains.  Future  tests 
with  the  segmented  silo  model  in  the  manner  described  here  should  obviously 
be  accompanied  by  more  flexible  soil  and  containment  vessel  wall  treatment 
to  reduce  the  shear  generated  at  the  soil  boundary  to  a  minimum. 
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SYMBOL  INDEX 


Section  Y 


a 


c 


crc  2 
F1’F2 


L 

n 


P 

P 

w 

W 

G 


P 

T 

0 


shell  radius 

cohesive  stress 

elemental  soil  cohesive  forces 

elemental  soil  friction  forces 

elemental  circumferential  soil  force  resultant 

silo  length  above  neutral  plane 

normal  stress 

surface  overpressure 

elemental  overpressure  force 

=  />L/2 

elemental  soil  wedge  weight 

failure  angle  through  soil  medium 

soil  density- 

shear  stress 

angle  of  internal  friction 
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VI. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  analysis  of  the  response  of  a  cylindrical  shell  to  a  plane  wave 
has  progressed  to  the  point  where  the  characteristics  of  the  predicted 
response  appear  to  be  in  excellent  qualitative  agreement  with  the  meager 
experimental  evidence  available  (e.  g.  NTS  tests  on  buried  conduits).  The 
results  obtained  with  a  step  pulse  loading  indicate  that  the  maximum 
stresses  and  displacements  may  be  significantly  lower  for  a  decaying  stress 
pulse.  This  conclusion  follows  from  the  fact  that  the  maximum  dynamic 
response  values  occur  in  the  neighborhood  of  five  or  six  transit  times.  Thus, 
it  is  necessary  to  investigate  effects  of  a  decaying  pressure  pulse  to  obtain 
realistic  results  for  full  scale  structures.  It  is  felt  that  a  logical  and 
important  continuation  of  this  phase  of  the  theoretical  investigation  should 
be  concerned  with  the  following  additions  or  modifications  to  the  existing 
method  of  analysis: 

1.  A  refined  stress  wave  with  a  linear  rise  to  a  peak  stress  and 
then  an  exponential  or  linear  decay  (depending  on  the  complexity  of  the 
Fourier  series  representation)  to  a  constant  stress  level. 

Z.  A  nonlinear  foundation  modulus. 

3.  Tangential  surface  shears  proportional  to  the  normal  pressure. 

In  view  of  the  fact  that  these  features  are  associated  with  real  soil  behavior 
the  possibility  of  including  these  items  should  be  thoroughly  investigated. 

The  ultimate  objective  of  the  experimental  investigations  on 
buried  cylindrical  shells  is  to  verify  the  ability  of  the  method  of  analysis 
to  predict  shell  stresses  and  displacements.  The  static  testing  of  tunnel 
models  represents  an  initial  step  toward  this  objective.  Unfortunately  the 
test  results  suffer  from  the  current  lack  of  satisfactory  soil  stress  and 
strain  measurement  devices.  Furthermore,  the  measurement  of  applied 
surface  pressure  and  related  surface  deflections  are  insufficient  (for 
the  pressure  vessel  used  for  these  tests)  to  accurately  determine  the 
required  free -field  stress  and  strain  at  the  tunnel  model  instrumentation 
section.  The  limited  experimental  results  of  Section  III  appear  to 
substantiate  the  static  theory  given  in  Section  II.  Two  hypotheses  are  re¬ 
quired  to  reach  this  conclusion;  namely,  that  1)  the  ratio  of  vertical 
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free -field  pressure  to  surface  pressure  decreases,  and  2)  the  ratio  of 
free -field  strain  to  average  strain  (based  on  the  surface  displacements) 
increases  with  increasing  containment  vessel  wall  friction.  It  is  recom¬ 
mended  that  these  tests  be  continued  as  soon  as  suitable  soil  pressure  and 
strain  gages  are  available. 

The  reported  theoretical  investigation  of  longitudinal  bending  of 
silo  linings  from  differential  horizontal  displacements  of  the  free -field 
relates  only  to  static  behavior.  That  is,  expressions  are  presented  to 
determine  lining  pressures,  shears  and  bending  moments  in  terms  of  a 
given  free-field  displacement  profile  and  foundation  modulus.  The  elastic 
stresses  produced  from  this  effect  are  shown  to  be  significant  for  the 
reinforced  concrete  silo  lining  studied.  This  investigation  should  be 
expanded  to  investigate  the  dynamic  effects  of  this  loading  source. 

For  such  an  analysis  it  would  be  necessary  to  obtain  a  time 
dependent  free-field  displacement  function  which  varies  with  depth.  The 
final  horizontal  displacements  could  be  made  to  conform  closely  to  those 
specified  in  the  design  criteria  for  current  ICBM  facilities.  It  is  anticipated 
that  the  horizontal  displacements  generated  from  Rayleigh  waves  would  be 
useful  in  specifying  the  displacements.  Two  structure -soil  interaction 
restoring  forces  should  be  considered,  one  proportional  to  the  relative 
displacement  and  the  other  proportional  to  the  relative  velocity.  Virtual 
mass  effects  should  also  be  included  in  the  analysis.  The  static  solution 
results  indicate  that  the  bending  rigidity  of  silos  is  so  great  that  only 
deflections  due  to  shear  need  be  taken  into  account.  This  suggests  the  use 
of  a  simplified  structural  model  for  the  investigation  of  dynamic  effects. 

The  exploratory  series  of  tests  performed,  to  investigate  the 
problem  of  longitudinal  silo  friction,  generated  by  relative  vertical  motion 
between  the  silo  and  the  soil  indicate  the  need  for  a  more  detailed  analysis 
along  the  following  lines.  As  a  first  approximation  to  the  analysis,  the 
zones  above  and  below  the  neutral  plane  can  be  considered  separately. 

Below  the  neutral  plane,  near  the  base  of  the  silo,  the  state  of  stress  is 
similar  to  that  of  a  deep  circular  footing,  for  which  the  shear  pattern  and 
corresponding  stresses  are  already  known  (13).  Above  the  neutral  plane, 
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solutions  can  be  obtained  for  the  plastic  stress  distribution  in  this 
axisymmetric  problem,  which  is  analogous  to  the  familiar  two-dimensional 
earth  pressure  theories.  Similar  solutions  have  been  developed  (14) 
previously  which  should  now  be  extended  for  this  problem.  These  solutions 
can  be  superimposed  to  yield  a  combined  solution  for  the  silo  problem. 

This  first  approximation  can  be  improved  by  considering  the  interaction 
of  these  solutions  to  satisfy  continuity  conditions  required  in  the  actual 
problem. 

Silo  model  push-out  tests  should  be  continued  employing  models 
with  various  length-to-r adius  ratios  and  alternate  soils  in  addition  to 
Ottawa  sand.  Tests  with  dynamically  applied  surface  pressure  should  also  be 
conducted  employing  a  segmented  silo  model  with  higher  axial  rigidity. 

The  technique  presented  for  the  analysis  of  the  stress -wave 
transmission  and  spallation  problem  in  a  silo  resulting  from  surface  air- 
blast  loading  indicates  that  severe  tensile  forces  do  not  develop  under 
loading  of  sufficiently  long  duration  compared  to  the  transit  time  of  the 
stress  wave  down  the  silo  length.  Further  investigation  should  be 
directed  toward  establishing  the  load  parameters  under  which  spallation 
damage  does  occur. 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


-  139  - 


REFERENCES 


1.  Newmark,  W.  M.  and  Hall,  W.  J. ,  Preliminary  Design  Methods 

for  Underground  Protective  Structures,  AFSWC-TR -60-5,  University 
of  Illinois,  Urbana,  Ill.,  December  1959. 

2.  Salmon,  M.  A.,  Investigation  of  Silo  Linings,  Final  Report,  Contract 
AF29(  60  1)  -  1 168  for  AFSWC,  Kirtland  AFB,  New  Mexico, 

Armour  Research  Foundation,  Chicago,  Illinois,  September  1959. 
Confidential. 

3.  Watkins,  R.  K.  and  Spangler,  M.  G.  ,  "Some  Characteristics  of  the 
Modulus  of  Passive  Resistance  of  Soil:  A  Study  of  Similtude",  Highway 
Research  Board  Proceedings,  37  (  1958)  . 

4.  Watkins,  R.  K.  ,  "Influence  of  Soil  Characteristics  on  Deformation  of 
Embedded  Flexible  Pipe  Culverts",  Highway  Research  Board  Bulletin 
223  (  1959) . 

5.  Timoshenko,  S.  ,  Theory  of  Elastic  Stability,  McGraw-Hill  Book  Co. , 

Inc.  ,  (  1936)  . 

6.  Bleich,  H.  H.  and  DiMaggio,  F.  L.  ,  "Dynamic  Buckling  of  Submerged 
Plates  and  Shells",  Columbia  University  Publication,  Dept,  of  Civil 
Engineering  Mechanics,  Office  of  Naval  Research  Contract  NONR-266(  08)  . 
Technical  Report  No.  12  September,  1954. 

7.  Salmon,  M.  A.  ,  Investigation  of  Silo  Linings,  Interim  Report, 

Contract  No.  AF29{  601) -2596  for  AFSWC,  Kirtland  AFB,  New  Mexico 
Armour  Research  Foundation,  Chicago/  Illinois,  October  I960, 
Confidential. 

8.  Wiedermann,  A.  H.  ,  The  Investigation  of  a  Concept  for  Medium- Structure 
Interaction,  Contract  AF29(  60 1) -2652  for  AFSWC,  Kirtland  AFB,  New 
Mexico,  Armour  Research  Foundation,  Chicago,  Illinois,  April  1961. 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


-  140  - 


) 


I 

I 

I 


9.  Anderson,  D.  C.  and  Ashby,  A.  ,  Development  and  Evaluation  of 
Instrument  Techniques  to  Measure  Structural  Effects  in  Nuclear 
Weapons  Effects  Sinulation,  Contract  AF  33(  616) -7098,  ASD 
Technical  Report  6 1-24,  Wright-Patterson  AFB,  Ohio,  Armour 
Research  Foundation,  Chicago,  Illinois,  April  1961. 

10.  Haywood,  J.  H.  and  Wilson,  L.  B.  ,  "The  Strain- Energy  Expression 
for  Thin  Elastic  Shells",  Trans.  ASME,  December  1958. 

11.  Spangler,  M.  G.  ,  "The  Structural  Design  of  Flexible  Pipe  Culverts ", 
Iowa  State  College,  Eng.  Exper.  Sta.  Bull.  153  (  1941)  . 

12.  Riley,  W.  F.  ,  Daniel,  I.  M.  and  Durelli,  A.  J.  ,  Stress  Wave 
Phenomena  in  Semi-Solids,  Contract  No.  AF29(  601) -2564  for 
AFSWC,  Kirtland  AFB,  New  Mexico,  Armour  Research  Foundation, 
Chicago,  Illinois,  July  1 9 6 1 . 

13.  Meyerhof,  G.  G. ,  "The  Ultimate  Bearing  Capacity  of  Foundations", 
Geotechnique,  2,  301  (  1951). 

14.  Sokolovski,  V.  V.  ,  Statics  of  Soil  Media,  Butterworths  Scientific 
Publications,  London,  (  I960)  . 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


141  - 


APPENDIX  A 

AIR-BLAST  INDUCED  STRESS  WAVES  IN 
LAUNCH  SILOS 

by  T.  A.  Zaker 


APPENDIX  A 


AIR-BLAST  INDUCED  STRESS  WAVES  IN 
LAUNCH  SILOS 

by  T.  A.  Zaker 

In  the  investigation  of  the  vulnerability  of  hardened  missile  launch 
facilities  to  nuclear  attack,  there  arises  the  possibility  of  damage  to  the 
underground  launch  silo  structure  resulting  from  transient  longitudinal  stress 
conditions  induced  by  air  blast  at  the  ground  surface.  For  certain  conditions 
of  loading  and  silo  wall  geometry,  this  phenomenon  can  result  in  the  develop¬ 
ment  of  severe  cracking  or  spallation  of  the  wall  structure.  The  purpose 
of  this  appendix  is  to  describe  a  brief  investigation  of  this  effect  in  a 
representative  launch  silo. 

In  this  appendix  a  technique  is  presented  for  an  analysis  of 
longitudinal  stress-wave  propagation  in  the  reinforced  concrete  walls  of  a 
missile  launch  silo  caused  by  surface  air-blast  pressure.  The  equations 
governing  the  phenomenon  are  given;  the  pertinent  physical  properties  of 
the  materials  are  discussed;  and  a  preliminary  analysis  of  an  idealized 
wall  configuration  is  presented  for  a  surface  loading  having  a  peak  over¬ 
pressure  of  300  psi  generated  by  the  surface  burst  of  a  1-MT  weapon. 

1.  Theoretical  Considerations 

When  a  column  of  material  is  subjected  to  time -dependent  surface 
pressure  loading  at  one  end,  longitudinal  waves  of  stress  are  induced  and 
propagate  down  the  length  of  the  column.  Such  waves  give  rise  to  strains 
and  particle  velocities  which  depend  on  the  properties  of  the  medium.  When 
the  stresses  exceed  the  proportional  limit  of  the  material,  the  waves  are 
propagated  at  velocities  which  depend  on  the  stress  amplitude  and  on  the 
current  physical  state  of  the  material.  Failure  by  spallation  occurs  when 
the  local  tensile  strength  c£  the  medium  is  exceeded,  and  the  material 
separates  locally.  The  problem  can  be  analyzed  on  the  basis  of  a  one - 
dimensional  model  of  the  system,  wherein  the  physical  state  variables  are 
regarded  as  functions  of  distance  along  the  column  and  of  the  time. 
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Consider  a  column  of  material  of  constant  cross  section,  and  assume 
that  the  stress  and  particle  velocity  components  parallel  to  the  column  axis 
are  uniformly  distributed  across  each  section  at  every  instant  of  time.  The 
differential  equations  for  the  local  conservation  of  mass  and  momentum 
are 


b  U 

bx 


C)U  b  CT 

Jo  dt  dx  ~ 


(A-  1) 


(A-2) 


where 

x  is  initial  position  coordinate  of  particle, 
t  is  time, 
p  is  density, 

is  initial  density,  a  constant, 
u  is  particle  velocity, 

o-  is  compressive  longitudinal  stress,  and 
f  is  body  force  per  unit  mass. 

It  is  to  be  emphasized  that  these  equations  are  written  with  respect  to  a 
Lagrangian  frame  of  reference,  in  which  the  initial  particle  position 
coordinate  (x)  and  time  (t)  are  taken  as  independent  variables;  no  limitations 
are  placed  on  the  magnitudes  of  particle  displacement  and  velocity.  This 
representation  is  particularly  convenient  when  the  current  physical  state 
of  a  particle  depends  on  its  prior  deformation  history.  If  we  supplement 
equations  A-l  and  A-2  with  constitutive  equations  of  the  medium  in  the  form 
of  stress-density  relations  [a  =  cr( )]]  ,  a  system  of  two  partial  differential 
equations  in  the  two  dependent  variables,  u  and  <r,  results. 

It  is  convenient  to  define  the  compressive  longitudinal  strain 
(*)  as  the  current  change  in  length  per  unit  of  initial  length.  Then  €  =  1  -  PQlp> 
and  if  we  assume  time -independent  stress -strain  relations  of  the  form 
it  =  <7  («),  equations  A-l  and  A-2  become 
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c)  u  b  o' 

b  t  b  x 
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where  E  is  defined  by  E  =  da  /de  and  is  not  necessarily  a  constant. 
Equations  A- 3  and  A-4  are  of  hyperbolic  type  and  can  be  written  in 
differential  form  along  so-called  characteristic  lines  in  the  x-t  plane  as 
follows: 


da  +  p  c  du  =  p  f  dx  (A-5) 

where  the  positive  and  negative  signs  apply  respectively  on  the  wave  lines 

dx  =  4-  c  d  t  (A-6) 

in  which  the  wave  speed  (c)  is  defined  by  c  =  J  E  / .  Strictly  speaking, 
c  is  simply  the  rate  at  which  a  wave  traverses  the  material  of  the  column; 
the  true  speed  of  a  wave  relative  to  a  frame  of  reference  fixed  in  space 
depends  not  only  on  c  but  also  on  the  local  particle  velocity  and  the  local 
strain. 

If  we  suppose  the  column  to  be  oriented  vertically  and  that  the  body 
force  (f)  is  due  to  gravity,  then  f  =  g,  where  g  is  the  gravitational 
acceleration.  We  may  then  define  a  =  a  -  p  g  x,  so  that  equations  A-5 
become 

da  +_  p  cdu  =  0  (A-7) 

on 

dx  =  +  c  d  t.  (A-6) 

These  are  the  usual  equations  for  the  one -dimensional  propagation  of  longi¬ 
tudinal  waves,  but  the  relations  between  a  and  €  here  are  functions  of 
position  in  view  of  the  definition  of  a  .  However,  the  body  force  resulting 
from  gravity  is  usually  assumed  small  compared  to  the  stress  gradients 
in  problems  of  interest. 
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Equations  A-6  and  A-7  lend  themselves  to  convenient  numerical  or 
graphical  computation  of  wave  motion  in  linear  or  nonlinear  media.  By 
considering  piecewise  linear  stress -strain  relations  and  piecewise  constant 
boundary  values  of  stress,  the  wave  motion  may  be  analyzed  on  the 
assumption  of  regions  of  uniform  physical  state  in  the  x-t  plane  separated 
by  steep-fronted  compression  or  unloading  waves  and  by  composition 
discontinuities.  The  analysis  of  each  wave -wave  interaction  (or  wave- 
composition  discontinuity  interaction)  reduces  to  the  application  of  the 
appropriate  one  of  equation  A-7  in  finite -difference  form  across  each  wave 
of  the  system  resulting  after  interaction,  and  the  satisfaction  of  the  interface 
conditions  of  continuity  of  stress  and  particle  velocity  in  zones  between  the 
waves.  Where  spallation  occurs,  the  condition  of  continuity  of  particle  velocity 
is  replaced  by  a  stress  boundary  condition. 

The  analysis  is  readily  extended  to  the  case  of  a  column  with  abrupt 
changes  in  cross  section  area*.  Retaining  the  one -dimensional  formulation 
and  neglecting  any  stress  transmitted  to  the  surrounding  medium  at  the  area 
changes  between  zones  of  constant  area,  the  condition  of  stress  continuity 
must  be  replaced  by  continuity  of  total  force  across  such  area  changes. 

Neglecting  body  forces,  equations  A-7  and  A-6  may  be  written  for 
each  constant-area  zone  as 

dF  +  pQc d  A  du  =  0  (A -8) 

on 

dx  =  +  c  d  t  (A-6) 


*  L.  Griffis,  The  Behavior  of  Longitudinal  Stress  Waves  Near  Discontinuities 
in  Bars  of  Plastic  Material,  NDRC  Report  No.  A-212,  OSRD  No.  1799, 
(September  1943). 
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where  F  =  <rA  is  total  compressive  normal  force,  and 
A  is  cross  sectional  area 

The  weighted  acoustic  impedance  (  c  A)  represents  the  local 
slope  of  characteristic  lines  in  a  plane  of  dependent  variables,  F  and  u. 

In  the  general  case,  and  A  are  functions  of  position,  and  the  local 

value  of  c  will  depend  on  the  current  state  of  the  material  if  the  medium 
is  nonlinear.  In  a  stepwise  numerical  computation,  these  values  are 
treated  as  locally  constant,  and  the  state  of  the  material  between  two 
opposite -facing  waves  resulting  from  a  wave -wave  interaction  or  a  wave- 
area  discontinuity  interaction  is  determined  by  the  intersection  of  two 
characteristics  in  the  F-u  plane.  These  characteristics  corrrespond 
respectively  to  the  two  opposite -facing  waves  arising  after  the  interaction. 

The  intersection  of  characteristics  in  the  F-u  plane  enforces  the  simultaneous 
equality  of  force  and  particle  velocity  in  adjacent  states  behind  opposite - 
facing  waves. 

To  illustrate  this  point,  consider  the  interaction  of  a  rightward- 
moving  wave  (dx  =  c^  dt)  in  zone  1  with  an  area  discontinuity  separating 
zones  1  and  2.  In  general,  the  interaction  will  give  rise  to  a  rightward- 
moving  transmitted  wave  (dx  =  c^dt)  in  zone  2  and  a  leftward -moving 
reflected  wave  (dx  =  -  c^dt)  in  zone  1.  Regions  of  uniform  state  in  the 
x-t  space  separated  by  the  discontinuity  lines  are  sketched  in  figure  A- la. 

The  first  digit  of  the  numerical  identification  denotes  the  constant-area  zone 
while  the  second  refers  to  the  physical  (F,  u)  state  shown  in  figure  A-lb. 
States  0  and  1  are  assumed  known,  and  are  connected  by  equation  A -8  with 
the  lower  sign  chosen  so  as  to  apply  across  the  rightward-moving  incident 
wave : 

FrF0  =  fo  ciAi  (uruo)-  {A_8a) 

Similarly  equation  A -8  with  lower  sign  applies  to  the  unknown  jump  in 
properties  across  the  rightward-moving  transmitted  wave,  while  equation 
A-8  with  upper  sign  applies  to  the  unknown  jump  in  properties  across  the 
leftward-moving  reflected  wave: 
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(A-8b) 


F,  -  F  =  P  c,  A,(u.  -  u  ) 
Z  o  o  Z  Z  Z  o 


F2  -  F1  --/°o  c1A1!“2  '  “l> 


(A-8c) 


State  2  is  thus  determined  in  F-u  space  by  the  intersection  of  oppositely 
directed  characteristic  lines  from  each  constant-area  zone,  corresponding 
to  equation  A-8  applied  in  finite  form  across  the  opposite -facing  resultant 
waves.  The  continuity  of  force  and  velocity  across  the  area  change  are 
obviously  satisfied  by  the  intersection  of  these  lines.  In  the  example 
sketched  in  figure  A-l,  the  transmitted  wave  is  a  wave  of  compression, 
while  the  reflected  wave  is  determined  to  be  an  unloading  or  expansion  wave. 

By  these  procedures,  the  solution  to  a  problem  of  one -dimensional 
stress  wave  propagation  can  be  carried  forward  graphically  in  tirrte  and 
distance  from  given  initial  and  boundary  value  data.  Specification  of  tensile 
fracture  criteria  leads  to  an  assessment  of  the  possibility  of  spallation 
in  the  material;  quite  complicated  failure  criteria,  such  as  time -dependent 
or  displacement-dependent  forms,  can  be  handled  by  a  semigraphical 
method. 


2,  Material  Properties 

The  cylinder  outside  walls  of  a  representative  launch  silo  are  of 
reinforced  concrete  extending  to  a  depth  of  about  163.  5  ft  below  grade,  and 
have  an  inside  radius  of  27.  5  ft.  The  wall  thickness  to  level  2  at  a  depth 
of  about  30  ft  is  8  ft,  and  the  remainder  of  the  wall  is  4  ft  in  thickness.  An 
annular  shoulder  of  inside  radius  13.  25  ft  and  approximately  4  ft  in  depth 
connects  to  the  silo  wall  at  grade,  and  a  circular  slab  7-ft  thick,  carrying 
the  flame  deflector  structure,  forms  the  base  of  the  silo. 

The  structural  concrete  used  is  of  3,  750-psi  ultimate  static 
compressive  strength.  Longitudinal  (vertical)  reinforcing  steel  is  specified 
as  No.  14  bars  spaced  at  9  inches  in  each  face  in  the  8-ft  thick  wall 
section,  and  No.  14  bars  spaced  at  18  inches  in  each  face  in  the  4-ft 
section.  Deformed  billet  steel  bars  are  called  for  in  accordance  with  ASTM 
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specification  A432-59T*.  These  bars  have  a  static  yield  strength  of 
60,  000  psi  and  an  ultimate  tensile  strength  of  90,  000  psi  at  7  percent 
elongation.  Each  bar  has  a  cross  sectional  area  of  2.  25  sq  in. 

Adynamic  yield  stress  of  4,875  psi  is  specified  in  the  concrete, 
and  a  minimum  dynamic  yield  stress  of  75,  000  psi  is  called  for  in  the 
reinforcing  steel. 

The  area  of  the  annular  shoulder  presented  to  surface  pressure 
(loaded  area)  is  3410  sq  ft.  The  silo  closure  is  carried  on  a  girder 
structure,  and  is  assumed  to  act  independently  of  the  silo  walls.  The 
cross  sectional  area  of  the  8-ft  section  (Section  1)  is  1582  sq  ft,  while  that 
of  the  4-ft  section  (Section  2)  is  742  sq  ft.  The  area  of  the  base  is  3120 
sq  ft. 

The  weight  density  of  concrete  is  taken  to  be  150  lb  per  cu  ft, 
while  that  of  the  soil  below  the  silo  is  taken  as  120  lb  per  cu  ft.  The  seismic 
velocity  in  the  soil  below  is  assumed  to  be  7000  fps. 

From  these  data,  composite  force-strain  curves  were  constructed 
for  sections  1  and  2  of  the  silo  by  addition  of  piecewise  linear  force-strain 
curves  for  steel  and  concrete  at  equal  strain.  The  total  steel  area  in 
section  1  is  1187  sq  in.  ,  and  that  in  section  2  is  556  sq  in.  assuming  6  inches 
of  concrete  cover  over  the  steel  circles  in  each  face.  Young's  modulus 
for  steel  was  taken  as  30  x  10  psi,  while  that  for  concrete  was  assumed 
to  be  1000  times  the  static  compressive  strength,  or  3.  75  x  10^  psi.  Perfectly 
plastic  behavior  was  assumed  in  the  concrete  at  compressive  strains  greater 
than  the  yield  strain  of  0.  001,  and  the  concrete  was  assumed  ineffective  in 
tension. 

The  piecewise  linear  force-strain  curves  for  concrete  and  steel 
in  sections  1  and  2  are  shown  in  figure  A-2,  while  the  composite  curves 
are  given  in  figure  A -3.  It  is  apparent  that  the  contribution  of  the  steel  to 
the  strength  of  the  composite  in  compression  is  quite  small.  On  the  other 

*  Deformed  Billet  Steel  Bars  for  Concrete  Reinforcement  with  60,  000  psi 
Minimum  Yield  Point,  ASTM  Stds.  ,  1959  Suppl.  ,  Pt.  1,  242-246. 
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hand,  the  strength  of  the  steel  in  tension  is  of  the  order  of  the  tensile 
strength  of  the  concrete,  if  the  latter  were  taken  at  the  commonly  assumed 
value  of  10  percent  of  the  compressive  strength. 

The  curves  shown  in  figure  A-3  can  be  applied  as  the  constitutive 
data  in  a  dynamic  problem  if  no  credit  is  taken  for  any  increase  in  strength 
owing  to  dynamic  effects;  such  increases  are,  at  best,  uncertain  quantities, 
and  their  omission  will  generally  result  in  an  underestimate  of  the  applied 
loading  required  to  cause  significant  damage. 

Hysteresis  effects,  resulting  from  compression  beyond  the  yield 
point  and  subsequent  unloading,  are  shown  in  the  individual  force -strain 
curves  and  in  the  composite  curves.  Hysteresis  effects,  in  general,  give 
rise  to  permanent  deformation  or  residual  strains  under  sufficiently 
intense  loading.  The  wave  speed  corresponding  to  each  linear  segment  of 
the  composite  force -strain  relations  is  to  be  calculated  from 


where  A  is  the  nominal  cross  sectional  area.  With  the  assumption  of 
negligible  tensile  strength  in  concrete,  it  is  seen  from  figure  A-3  that 
severe  yielding  and  steel  fracture  would  occur  at  tensile  forces  of  about 
110  x  10^  lb  in  section  1  and  about  50  x  10^  lb  in  section  2. 

3.  Stress  Wave  Problem 

In  view  of  the  discussion  of  material  properties,  some  aspects 
of  the  stress  wave  propagation  problem  in  silo  walls  can  be  illustrated  by 
investigation  of  a  simplified  model  neglecting  the  effect  of  reinforcing 
steel.  To  this  end,  we  consider  a  model  of  the  silo  wall  consisting  of  a 
column  of  plain  concrete  having  a  single  area  change  (Sections  1  and  2), 
and  assume  the  tensile  strength  of  concrete  to  be  10  percent  of  its  static 
compressive  strength,  or  375  psi. 
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The  masses  of  the  annular  shoulder  and  base  slab  are  neglected,  and 
these  components  are  assumed  to  act  simply  as  load  transmitters  from  the 
air  and  to  the  soil  below.  The  force  at  the  base  of  the  silo  walls  is  assumed 
to  be  transmitted  to  a  semi-infinite  column  of  soil  (Section  3)  equal  in  area 
to  the  silo  base  slab,  and  linear  elastic  behavior  is  assumed  in  the  soil. 
Properties  of  sections  1,  2,  and  3  for  the  elastic  range  of  behavior  are 
given  in  the  table  A-l.  The  wave  speed  (c)  for  the  concrete  is  calculated 

from  c  =  y  E I pQ. 


T able  A-l 


Section 

0 

i 

2 

3 

Component 

- - - L 

Shoulde  r 

_ 

8 -ft  wall 

4 -ft  wall 

Soil 

- L 

A,  sq  ft 

1 

1 

3410 

_ 

— 

1582 

742 

3120 

2  4 

p  ,  lb- sec  /in.  1 
To 

_ 

1 

| 

J _ 

2.  25xl0-4 

2.  25xl0"4 

1.  80x10  4 

IB 

10.  65xl03 

10.  65xl03 

7.  OxlO3 

p  cA,  lb-  sec  /in. 
Jo 

— 

6.  62xl06 

3.  10xl06 

6.  78xl06 

The  loading  is  assumed  to  be  a  pulse  of  peak  overpressure  300  psi 
resulting  from  the  surface  burst  of  a  1-MT  weapon.  A  simplified  representa¬ 
tion  of  the  overpressure -time  history  at  the  ground  for  surface  busts  has 
been  developed  by  Newmark*  from  theoretical  solutions  to  the  point- source 
explosion  problem  in  free  air.  The  representation  is  based  on  the  replace¬ 
ment  of  the  actual  pressure-time  curve  by  two  straight  line  segments.  From  the 


*  N.  M.  Newmark  and  W.  J.  Hall,  Preliminary  Design  Methods  for  Under¬ 
ground  Protective  Structures  (U),  AFSWC  TR-60-5  p.  3-16  (December, 
1959)  SECRET. 
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initial  peak  of  300  psi,  the  overpressure  decays  linearly  with  time  to  about 
68  psi  at  time  t^  =  62  msec  and  then  linearly  to  zero  at  time 

t^  =  700  msec. 

In  view  of  the  area  of  the  shoulder,  the  initial  peak  force  applied 
to  the  system  is  147  x  10^  lb.  The  wave  pattern  calculated  by  the  techniques 
presented  earlier  is  shown  m  the  x-t  diagram  of  figure  A-4  up  to  t  =  35  msec, 
together  with  a  schematic  diagram  of  the  model  sketched  along  the  distance 
axis.  The  linear  decay  of  the  loading  is  replaced  as  shown  by  two  constant- 
pressure  intervals,  the  decrease  from  the  initial  pressure  occurring  at 
about  t  =  11.  5  msec.  The  initial  digit  of  the  numerical  identification  of 
regions  in  the  x-t  diagram  of  figure  A-4  refer  to  the  constant-area  zones 
in  the  model,  while  the  remaining  digits  refer  to  the  corresponding  force - 
velocity  states  shown  for  clarity  in  two  F-u  diagrams,  figures  A-5  and 
A-6.  Compression  waves  are  shown  in  figure  A-4  as  solid  lines,  while 
unloading  waves  are  shown  as  dashed  lines. 

The  principal  characteristics  through  the  origin  of  the  F-u  plane 
for  the  linear  elastic  range  in  each  section  of  the  concrete  column  and  in  the 
soil  are  shown  in  figures  A-5  and  A-6.  The  slope  of  each  of  these  lines  is 
equal  to  the  corresponding  weighted  acoustic  impedance  (pqcA).  Every  other 
line  segment  in  these  diagrams  has  a  slope  equal  in  absolute  value  to  that 
of  the  corresponding  principal  characteristic  through  the  origin.  The 
tensile  forces  required  to  produce  failure  in  sections  1  and  2  are  also  shown 
along  the  negative  F-axis  on  the  assumption  that  the  tensile  strength  of 
concrete  is  10  percent  of  its  static  compressive  strength. 

It  is  apparent  that,  for  the  surface  loading  considered  in  this 
problem,  the  material  remains  elastic  and  well  within  acceptable  stress 
limits.  Velocities  up  to  about  40  ips  are  attained  (State  4,  Fig.  A-5),  and 
a  maximum  compressive  force  of  about  180  x  10^  lb  develops  (State  16, 

Fig.  A-5)  at  the  base  of  the  silo  wall  upon  reflection  of  waves  at  the 
silo-soil  interface.  As  the  surface  pressure  loading  continues  to  decay, 
the  forces  and  velocities  in  the  system  can  be  expected  to  decrease 
continually  toward  zero  (Fig.  A-6). 

From  figure  A-4,  observe  that  the  surface  pressure  decays  to 

half  its  initial  value  at  about  t  =  40  msec.  This  characteristic  of  the  loading 
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corresponds  to  nearly  three  transits  of  a  stress  wave  down  the  length  of 
the  silo.  Owing  to  the  relatively  slow  decay  of  the  assumed  surface 
loading,  no  appreciable  tensile  stresses  are  expected  to  develop  in  the 
model.  On  the  other  hand,  it  is  clear  that  for  intense  surface  pressure 
of  sufficiently  short  duration  (of  the  order  of  a  fraction  of  the  transit  time 
of  a  stress  wave  down  the  silo  length)  severe  tensile  forces  may  develop  near 
the  area  discontinuity  between  sections  1  and  2  upon  emergence  at  the  surface 
of  the  reflection  of  the  inital  input  wave  from  the  base  of  the  silo  (note 
State  21,  Fig.  A-6).  This  indicates  that  spallation  damage  to  the  silo 
structure  is  more  likely  to  occur  as  a  result  of  intense,  short-duration 
loading  from,  say,  the  air  blast  caused  by  a  close  near-surface  high- 
explosive  burst  than  from  that  caused  by  MT -range  nuclear  weapon  in  the 
300-psi  overpressure  region.  It  is  precisely  this  effect  which  accounts  for 
the  so-called  close-in  "brisance"  of  a  high  explosive  and  leads  to  severe 
scabbing  and  fracturing  of  solids  in  contact  explosions. 

Finally,  it  may  be  noted  that  body  forces  have  tacitly  been  assumed 
negligible  in  the  foregoing  analysis.  The  dead-load  compressive  stress 
at  the  base  of  the  silo  wall  may,  in  reality,  be  as  high  as  about  200  psi  if 
the  entire  weight  of  the  silo  wall  structure  is  carried  by  the  walls  and  base 
slab.  However,  the  material  is  found  to  remain  linearly  elastic  throughout 
the  motion  considered  in  this  problem,  and,  in  such  a  case,  it  can  be  shown 
that  dead-load  forces  may  simply  be  added  to  those  obtained  from  the 
analysis  neglecting  body  forces,  in  accordance  with  a  transformation  similar 
to  that  discussed  in  connection  with  equation  A-7.  In  general,  this  is  not 
possible  when  plastic  flow  and  fracture  occur  in  a  problem. 

4.  Conclusion 

A  technique  has  been  presented  for  an  analysis  of  the  stress-wave 
transmission  and  spallation  problem  in  a  representative  missile  launch  silo 
resulting  from  surface  air-blast  pressure  loading.  A  preliminary  analysis 
of  the  problem  has  been  carried  out  on  the  basis  of  a  highly  simplified  model 
of  the  silo-wall  system,  and  composite  force -strain  curves  for  the 
reinforced  concrete  of  the  silo  walls  have  been  devised  for  use  in  a  more 
sophisticated  analysis. 
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The  analysis  of  the  plain  concrete  model  under  loading  resulting 
from  a  1-MT  surface  burst  in  the  300-psi  overpressure  region  indicates 
that  decay  times  associated  with  this  type  of  loading  are  sufficiently  large 
compared  to  the  transit  time  of  a  stress  wave  down  the  length  of  the  silo  so 
that  severe  tensile  forces  are  not  likely  to  develop.  On  the  other  hand,  it  is 
observed  that  intense,  short-duration  air-blast  pressure  may  result  in 
strong  tensile  waves  and  consequently  spallation  or  scabbing  of  the  material. 

Further  investigation  should  be  directed  toward  the  establishment 
of  load  parameters  under  which  severe  damage  of  this  type  may  occur.  In 
addition,  significant  improvement  can  be  made  to  the  one -dimensional  model 
described  here.  These  improvements  should  be  based  on  more  realistic 
treatment  of  coupling  between  the  motion  of  silo  walls  and  that  of  the  closure 
and  base  slab,  and  between  the  motion  of  the  base  slab  and  that  of  the  soil 
below  the  silo.  Various  degrees  of  elastic  and  anelastic  coupling  can  be 
considered  along  with  equivalent-mass  addition  to  the  system  with  time, 
within  the  frame -work  of  the  one -dimensional  formulation. 
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COMPOSITE  FORCE -ST RAIN  RELATIONS  FOR  SILO  WALLS 
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PLAIN  STRAIN  SOLUTION  FOR  THE  DISPLACEMENTS  AT 
BOUNDARY  OF  A  CIRCULAR  HOLE  IN  AN  INFINITE 

ELASTIC  PLATE 


The  plain  strain  solution  for  the  displacements  at  the  boundary  of 
a  circular  hole  in  an  infinite  elastic  plate  for  the  case  of  a  radial  load  varying 
as  cos  n  9  can  be  obtained  from  the  general  solution  of  the  two-dimensional 
problems  in  polar  coordinates*.  The  stress  function  for  loads  varying  as 
cos  n  9  is 


i  f  n  .  -n  2  +  n  ,  2  -n 

9  =  [_a-r  +  b„r  +  c_r  +  d„r 


n  n 


n 


sin  n  9 
cos  n  9 


(B  - 1) 


where  n  ^  2 

The  corresponding  stresses  are 


2.  n-2 


2.  -n  +  2 


rr  l.  n 


a^(n-n  )r“  “  -b^  (n  +  n")  r  +  c  (2  +  n-n^)  r 


2.  n 


n 


n 


+  d^  (2-n-n^)  r  cos  n  9,  (B-2) 

2  .  n-2 


99 


=  a^(n  _n)  r  +  b^(n  +  n  )  r  **  '  *'  +  Cn(2+3n  +  n")  r 


+  d^(2  -  3n  +  n^)  r  cos  n  9, 


(B-3) 


2.  n 


r9  =  [an^n  "  r  *  **  "bn^n  +  r  "  “  +  cn(  n  +  O  r 


d  (n  -  n^)  r  n  1  sin  n  9  . 
n  -1 


(B-4) 


*  Timoshenko  and  Goodier,  Theory  of  Elasticity.  2nd  Edition,  p.  116, 
McGraw-Hill  Book  Co.  ,  N.  Y.  (1951). 
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We  consider  the  infinite  plate  with  a  circular  hole  of  radius,  a,  subjected  to 
the  internal  pressure,  p  cos  n  0,  where  n  ^  2.  The  boundary  conditions  at 
r  =  a  are 

<r  =  -p  cos  n  0,  ct^q  =  0,  at  r  =  a  .  (B-5) 

Since  the  loading  is  self-equilibrated,  the  stresses  at  infinity  must  go  to 

zero  as  — .  Hence 
r 


a  =  c 
n  n 


(B-6) 


Substituting  the  expressions  for  the  stresses  into  the  boundary  conditions, 
we  obtain 

-bn(n  +  n^)  a  n  ^  +  d^  (2-n-n^)  a  n  =  -p 

-b  (n  +  n^)  a  n  ^  +  d  (n-n^)  a  n  =  0  . 
n  n 


Hence 

b. 


n 


n 


pan+Z 

2(n+l) 

n 

2(n-  1) 


(B  -7) 


The  radial  displacement  at  r  =  a  for  plain  strain  is 


u 


=  |bn(l  +  i))  n  a ’ n "  1  +  dn  [4(1-  -\)Z)  -  (1  +1))  (2-n)]  a'n+1j  cos  n 


0 

(B-8) 


or,  using  equation  B-7, 


u 


pa 


r  =  a 


2n+  i 

L  ”2^7 

n  1 


(1+-0  )  (1  -  2  i)  )  J  cosn0 


(B-9) 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


B-2 


Now  —  -j  ^  ^  j - 2'^  y  is  the  apparent  modulus  of  uniaxial  strain  which  we 

have  defined  as  E 
value  of  k,  hence 


have  defined  as  Ec>  The  ratio  of  radial  pressure  to  displacement  gives  the 


ka  =  E 


2  i 

n  -  1 
2'n+"l 


,  n  2. 


(B-  10) 
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APPENDIX  C 


COMPUTER  PROGRAM  FOR  SHELL  RESPONSE  TO  RADIAL  PRESSURE 


The  flow  diagram  for  the  program  machine  code  is  given  in 
figure  C- 1  .  An  algebraic  programming  language  is  employed  based  on 
the  Internal  Translator  (IT)5*  originally  developed  at  Carnegie  Institute 
of  Technology,  and  serves  as  basic  input  to  the  UNIVAC  H05  compiler. 

The  following  nomenclature  has  been  used  in  the  computer  program.  Var¬ 
iables  marked  with  a  single  asterisk  are  required  input  data  for  the  program 
and  those  with  a  double  asterisk  are  output  variables.  Variables  defined 
within  the  program  are  not  listed  here. 


##  ct 

YO  =  Z  ,  ( — ),  nondimensional  time 

cL 

an 

Y1  =  x  ,  ( — ),  nondimensional  modal  displacement 
n  a 


Y2 

X  'n’  H 

Y20  = 

ka/ f 

** 

Y21 

s  !  j>  c 

Y22 

K  , 

m 


ct 

*  **  1 

C41  ,  Y23  =  Tj, 

#  jjcjje 

C40  ,  Y24  =  jo,  free-field  stress  biaxial  ratio 

jj{  jjc  sjc 

C43  ,  Y25  =  last  value  of  Z  (  T_  4  99  AlJ 


❖  *  # 

C30  ,  Y31  =  a/h,  shell  radius  to  thickness  ratio 


*  A.  J.  Perlis,  J.  W.  Smith,  and  H.  R.  Van  Zoeran,  "Internal  Translator 
(IT)  A  Compiler  for  the  650",  Computation  Center  Report, 

January  1957. 
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P'P  s 

E  /ka 


Y  32 

Y  33 

Y  34 


** 

** 

** 


Hf  ** 

C42  ,  Y35  = 

Y  49  = 


Y  104 


(T 

o  a 

TT  IT 

Ar^,,  Type -Out  interval 

uj^,  natural  frequency  in  nth  mode 

=  ka/E* 


C44*,  Y 105**  =  XA 
4 


1,  for  addition  of  fictitious  buckling  pressure 
0,  for  omission  of  fictitious  buckling  pressure 


CO 


** 


C31 

* 

C32 

* 

C34 

♦ 

C35 

* 

C36 

♦ 

C37 

♦ 

C38 

C39 

* 

Ill 

$ 

119 


124 


=  Sr , 

= 

-  Ps> 

=  C, 

=  s  , 


=  ka, 
=  E, 


integration  interval 

2  -4 

soil  density  (lb-sec  -in  ) 

2  -4 

shell  density  (lb-sec  -in  ) 
stress  wave  velocity  (ips) 

_3 

damping  coefficient  (lb-sec-in  ) 
peak  free -field  pressure  (psi) 
apparent  soil  modulus  (psi) 
foundation  modulus  times  shell  radius  (psi) 
Young's  modulus  of  shell  (psi) 
initial  displacement  mode 
final  displacement  mode 


0,  Short  Program,  No  output  at  9  =  ^  .  .  . 

[l,  Long  Program,  xn  and  p^  output  at  9  =  .  .  . 
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MJ1  Switch  < 


Off,  Typeout:  Header,  Rigid  Body  Mode  (n  =  1) 

Static  Solution,  and  ^  output 

[On,  Typeout:  Same  as  per  off  position  plus  each  mode 
output 


Table  C-l  identifies  the  four  columns  of  output  from  the  computer 
program.  Rows  that  do  not  have  the  letter,  U,  in  the  comments  column  are 
not  typed  if  the  MJ1  switch  is  off.  Unless  indicated  to  the  contrary,  all 
response  output  is  for  the  head-on  point  (0=0°).  If  the  input  value  of 
IZ4  =  0,  the  displacements  and  pressures  at  S  =  TC/4,  TC /2  .  .  .  are  not 

r— i 

computed.  In  the  static  and  summation  (  ^  output)  portion  of  the  list, 
the  following  abbreviations  are  used 


n. 


W(0) 


-I 


N 


I  N 

cos  n  8,  'j 

!  N 


2  for  n.  =  1 

l 

n.  otherwise 
l 


nF 

p(6)  =£ 

N 


r 


N 


cos  n  8, 


0  for  n .  =  1 

n.  otherwise 
i 
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t 


V 


SOLVE  EQUATION  OF 
MOTION  AT  JTME  r 
(  EVALUATE  a  ,  a  ) 

n  n 


> 


EXIT  NUMERICAL 
INTEGRATION  SUL  ROUTINE 


Fig.  C-l  FLOW  DIAGRAM  FO'FTM'A'CHINFr  CODE 
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Fig.  C-l  (Cont.  )  FLOW  DIAGRAM  FOR  MACHINE  CODE 


Tabic  C-  1 
OUTPUT  FORMAT 


Field  1 


ka/  pc 


Field  2 


sty  c 


Field  3 


K 


Field  4 


Comments 


W(0) 

W(X/4) 

P(TC/4) 


W(-rc/4) 

P(Tt/4) 


i  v 

1  ka/E 


i  c  a/kah 
o 


0 

6 


M  /c  a 
n  o 


M  /cr  a 
n  o 


<W/,To:  M/aoa 

W(tr/2)  W(3TC/4) 

P(-ir/2)  P(3lC/4) 


W(TC/2) 

PCK/2) 


|  Mn/,Joa 
j  W(3TC/4) 
P(3'ir/4) 


C1 

u 

^  ! 

j 

u 

P'Ps 

u 

Att 

u 

hr 

o 

p  /o' 
o  o 

Pl'% 

u 

Po/a; 

Pl^o 

u 

o  'T* 

0 

0 

0 

0 

p  fa 

n  o 

0 

p  fa 
rn  o 

0 

u 

0 

u 

P(0) 

u 

W(rc) 

u 

P(ic) 

u 

0 

U 

0 

u 

0 

u 

P(0) 

u 

1 

W  (1C) 

u 

p(‘0 

u 

Header 


Repeat 


r  Repeat  to 

J  rF 


Repeat 

to 


For 
n  =  1 


For  niil 
Repeat 
to  n_ 


Static 

Solution 


z 

Output 
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IT  CORRECTION  PROGRAM  00 a  u?3 

CHARGE.  <199*  OPTIONAL  OYNjAmIC  RING  BFNQ I NG  t  RRK  *  21 1 8 
I  3o  Y  530  C  1000  S 


SE?  11  1961 
HESL  3,366 
U21 

95  E2 1 3 . 4  *  6 » 8 


1 . 

57 

50.120* 130. 1 .999. 

F 

1.05 

65 

0,50  IF  120*530 

P 

1.1 

Y  120 ; »0. 0 

P 

2. 

50 

C I 20  s hO'O 

2.5 

C20HM.  1415927 

r 

2.7 


10  ttt2 


3. 

G2l 

3.1 

1 

Y20 t «C 3B /C 3 1*C 34*034 

3.2 

Y2 1 t«035/C31*C34 

3.3 

Y22  s  »r 33 

3.45 

Y23:»r4l 

3.5 

Y24  sar4o 

3.55 

Y25S«043 

3.6 

Y31 :»r30 

3.65 

Y32 ; »03l /C32 

3.7 

Y335S039/C33 

3.75 

Y34  j«(-36‘<'C30/C33 

3.6 

Y35 • »C42 

3 .85 

Y104  t«C3b/C37 

3.9 

Y 105  t  tfC44 

4. 

91 

TY20  TY21  TY22  TY2j 

4.1 

79 

TY104  TY105 

5. 

80 

TY24  TY26  TY31  TY>2 

6 . 

81 

TY33  rY‘j4  TY35 

r. 

92 

Y36!»Y2l/Y20 

a. 

Y305»I11*IU 

9. 

Y4 1 J«Y30/Y30>1 .0 

10. 

Y42:«y30-1.0 

11. 

Y37 :«y41*Y32»Y34*Y20 

12. 

Y30:f»Y4l**>Y3l*Y32 

13. 

Y38 ;8y30*Y21 

14. 

Y43:»y30*Y20 

15. 

V39;«y42»Y33/l2.0*Y3l*Y3l*Y3: 

16. 

Y40  ;Fy106¥Y33'KY34''2.0*Y31*Y3: 

17. 

Y47:«i .0+Y30*Y22/X j 1 

ia. 

Y48:0<  1 ,0+Y42*Y39 > *Y43/Y47 

19, 

Y49!«06E» Y460 

20. 

Y30s»r20/6.0*Y49 

21. 

Y26:»Y3S 

22. 

I 185«1 

23. 

I  15!Fi .0+Y35/Y30 

24. 

C0i«Y-<5/ll5 

26. 

T ( 300 . ) TOO 

26. 

C10J«i  .0  +  Y24 

27. 

C  11  :»<C1 0/2.0  J  ■*■(  1  •0-’Y24  )  /4.0 

26. 

C12:»( 1.0-Y24 1/2.0 

29. 

C13t»f 1.0-Y241/4.0 

30. 

C14t»r20*Y23 

31. 

C15:«t .0+Y104+Y24 

32. 

C 16 ! #y23/2. 0 

33. 

C17s«( 1,0+Y104-Y24)/2.0 

34. 

C16?«r  11+0.375 

33. 

C  19 ;»( 3.0-2. 0*Y24  >/B,0 

36. 

Y27 :Fy31*Y33/2.0*Yj3+Y31 

e 


e 

F 


e 


e 

F 

F 

F 

F 

F 

F 

F 


F 

r 


F 


F 

F 

r 

F 

F 
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it  n.  ii 


yt. 

Y2O»0Y42*Y33/12.O*Y34*Y3l*Y3l 

7 

37.1 

G76  IP  12400 

P 

37.2 

C25*0J 11*C20/4.0 

r 

37.3 

66 *123*1* 1*3* 

7 

37.4 

66 

C ( 20  + 123  > J«Q4£. 123*0250 

7 

37.5 

G70  Ip  111*1 

F 

37.6 

C24  *0-1.0 

F 

37.7 

G76 

7 

37.6 

78 

024  s  0-C24 

7 

38. 

76 

Y79  s  P Y27/Y3 1 

F 

39. 

Y29:»a.0*Y28+Y31*Y31 

7 

40. 

I 17 •  0?5 

7 

40.5 

I22*0?4 

r 

41. 

23 

Y0  *00.0 

F 

42  • 

Y 1 :«0.0 

r 

43. 

Y2 :»0 .0 

7 

44. 

T  (0.  ) 

7 

45. 

TUI  TY49  T(  0.  )  T  (  0  •  ) 

7 

46. 

TYO  TYl  TY2  KO.) 

F 

47. 

G7 

7 

48. 

2 

G35  IP  Y23//0 •  0 

7 

49. 

G i  17  iF  YO'02.0 

7 

50. 

Y44  J«n8E* 1 .O-YOQ 

7 

51. 

G29 

r 

52. 

25 

Y44:»C20 

7 

53. 

I  17:0-*, 4 

F 

54. 

29 

Y46 • 0 V44 

F 

55. 

I 12 i«t  11 

7 

56. 

I l4*0n 

7 

57. 

33 

G30  If  11200 

7 

56. 

Y( 60+J14) J0Q3E* I 1 2*Y460/ 112 

7 

59. 

G31 

7 

60. 

30 

Y ( 60+ 114) S0Y46 

7 

61. 

31 

32*1 1a *1*1*3* 

7 

61.5 

27 

I 1 3  I 0aO+1 14  +  116 

r 

62. 

G26 I F ( 112-116)00 

r 

63. 

VI  13  « rt (Q3E.  (  11 2-1 16  >*Y46Q/I  12-116) 

64. 

+Q3E* ( I12+I 16)*Y46C/i 12+1 16 

7 

65. 

G32 

7 

66. 

26 

Yl 13 ! 0Y45+Q3E ’ ( l 1 2+i 1 6 ) *Y46Q/1 12+1 16 

7 

67. 

32 

1 1  :it 1 1 

7 

68. 

G38  IF  11405 

7 

69. 

G3C  IP  114010 

P 

70. 

G39  Ip  114015 

7 

71. 

I 14*0)0 

r 

72. 

112*00 

7 

73. 

Y46J0Y44 

p 

74. 

G33 

P 

75. 

34 

G35  IF  YO • Y23 

C 

76. 

Y69*0i .O-YO 

7 

77. 

Y80s«(-Y69*C10*Y70)+jCl l+Y71)-( Y69*Cl2*Y72) 

78. 

+C 13*y73 

F 

7S.  1 

Y30:«(-Y69*Y70)+<  0.5*Y71*0.75+Y69*Y69>-(0.5*Y69*Y72 ) 

73.2 

+0. l2p*Y73 

F 

79. 

Y84*0y7O-<  Y69*Y7l >+0.5*Y72 

7 

79.5 

Y80 : 0 ( Y80+Y104*Y30+Y36 yY84 ) /C 14 

F 

80. 

Y81  ;0<  (-769*0 10 *Y6o)+(Cil*Y6l)-<Y69*C12*Y62) 

81. 

+C 13*y63 ) /C 14 

7 

82. 

Y82  * 0 ( (-Y69-Y60 )  +  ( 0*5*1 0,75  +  Yb9*Y69) *Y6l ) - ( 0* 5*YA9*Y62 ) 

83. 

+0. 12f*Y63)/C14 

p 
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83.  1 

Y82S»Yl04*Ytl2 

84. 

Y83S8,  Y60-< Y69  +  Y6 l 1 +0 . 5*Yb2 ) *Y 104/C  14 

S3. 

036 

86. 

35 

G58IF?.0'*»Y0-Y23 

86.1 

Y45l«r20 

Sfc.2 

G59 

8b. 3 

58 

Y45  saoSE*  I.O-YO+Y23Q 

87. 

59 

Y46 S8y45 

88. 

I 1258)11 

89. 

11458s 

90. 

G33 

91. 

38 

I 12580 

92. 

1  14  58,5 

93. 

G33 

94. 

39 

Y69S8, .0-Y0 

98. 

G22  IF  Y23&0 

96. 

4  1  •  I  1  p  « 0  '  1 »  3 « 

9/. 

40 

Y ( 90+ 1 12 ) 5# ( Y ( 60+1 12 ) "Y ( 65+1  l2)  ) /Y23 

98. 

41 

Y(  100  +  112)  |8(Y<70+U2>-V(  /5  +  U2)  I/Y23 

99. 

G56 

100. 

22 

28* I 1?*0« 1 1 3 • 

101. 

37 

Y ( 90  +  T 1 2  >  5»0«0 

102. 

20 

Y ( 100  +  112  >  j«0.0 

103. 

56 

Y80S8y75-( Y76*Y69 +C 1 6 > + ( Y77*0*5> 

103.1 

-( Y69*Y100)+(0.5*Yi01*0.75+Y69*Y69> 

103.2 

Y80  5  8yB0- <  0*  5kY69’I'y1  02  1  +0  •  1  25*Y  103 

104. 

Y30  5  8  (  C10*Y75 )  +  ( C  l2>itY77  )  -  (Y69*C  10*YloO) 

104.1 

+(C11*Y101 )-< Y69*Ci2*Y10? )+Cl3+Yl03 

105. 

Y64  5  5(Y76+YiOO-(  Y69*Y101  )+0.5*Yi02 

105.5 

Y00S8( Y30+Y104*Y30+Y36*Ye4)/c2o 

10b. 

Y61  58,  <C10>?Y65  J  +(Cj2yY67)-(  Y69*C10»Y<?0) 

107. 

+  (CU*Y91 )-( Y69*cl2«Y9H>+(C13*Y93> >/C2o 

ioa. 

Y82  5  8 Y65- ( < Y69+C 1 6  >  +Y66 ) + ( 0 . S+Y&7  > - l Y69+Y90 ) 

109. 

+0.5*,0.75+Y69wY69)*Y91 

110. 

Y32  58, YQ2-(0.5*Y69wY92l+0* 125*Y93) /C20 

iio.l 

Y8258Y104+V02 

111. 

YS3  5  8 , Y66+Y90- ( Yb9,Y9 1 ) +0 . 5*Y92  > /C20 

111.5 

Y6358YB3--Y104 

112. 

36 

Cl  555 Y? 

1 1 2 . 1 

G93  IF  111*1 

1 12.2 

C2  s8Yql+v83*C35*e34/C3B 

112.3 

C2  58  (  r36**C 2  )  -C 35  *C’34*C  1 

112.4 

C2 : 8C  3 /( < 2. vC32/COo>  +C33  +  C31 )»C34+C3« 

112.5 

54 

112.9 

93 

105810 

113. 

C258<  ,  (  YUl+Y62+Ys3*Y36  i*Y37> -(  V38+Y2! 

114. 

- ( 1 . 0+Y42*Y39-Y40 *y60  >  wYu3*Y  1 ) /Y47 

115. 

54 

116. 

3 

G5  IF  <Y26-0.5*CO>*Yo 

117. 

53 

G63  IF  111*1 

118.  • 

C300  5  e Y27*YB0 

119. 

C5C0 5tfY79*Y80 

120. 

TY26  TC300  T(0.)  Tc500 

120.  1 

C  (300  +  118)  S55C300 

120.2 

C (500+118) i«C500 

120.3 

577  IF  12480 

120.34 

125  555600  +  1  18 

120.38 

12658<}00  +  i  18 

120.4 

68.123* 125.  100.  12  b . 

120.5 

68 

C  1 23  5  t/C5O0 

e 
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120.6 

77 

C500<4Y81+C35*C34*( < Y83/Y 104 ) ”Y2*C 37/C36)/C37 

F 

120.63 

94 

TY i  Ty2  KO.I  TCSOo 

r 

120.66 

G95 

f 

120.7 

67 

C500««Y31+Y02+<  Yfl3*Y36)-<  Y 1 +Y36*Y2 > *Y3l/Y34 

r 

120.8 

TY 1  TC300  TC400  TC5OO 

F 

120.85 

95 

I0:4lo 

z 

120.9 

C ( 500+1 1 3 ) :4C ( 500+118 )  +C500 

r 

121. 

064  IP  11141 

F 

121.1 

C ( 400  +  1 1 8 ) :4C<400+I18)+C«00 
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